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ACTIVE  SHIELDING  AND  CONTROL  OF  ENVIRONMENTAL  NOISE* 

J.  LONCARICt,  V.  S.  RYABEN’KII^,  AND  S.  V.  TSYNKOV§ 


Abstract.  We  present  a  mathematical  framework  for  the  active  control  of  time-harmonic  acoustic  dis¬ 
turbances.  Unlike  many  existing  methodologies,  our  approach  provides  for  the  exact  volumetric  cancellation 
of  unwanted  noise  in  a  given  predetermined  region  of  space  while  leaving  unaltered  those  components  of 
the  total  acoustic  field  that  are  deemed  as  friendly.  Our  key  finding  is  that  for  eliminating  the  unwanted 
component  of  the  acoustic  field  in  a  given  area,  one  needs  to  know  relatively  little;  in  particular,  neither 
the  locations  nor  structure  nor  strength  of  the  exterior  noise  sources  needs  to  be  known.  Likewise,  there 
is  no  need  to  know  the  volumetric  properties  of  the  supporting  medium  across  which  the  acoustic  signals 
propagate,  except,  maybe,  in  the  narrow  area  of  space  near  the  boundary  (perimeter)  of  the  domain  to  be 
shielded.  The  controls  are  built  based  solely  on  the  measurements  performed  on  the  perimeter  of  the  region 
to  be  shielded;  moreover,  the  controls  themselves  (i.e.,  additional  sources)  are  concentrated  also  only  near 
this  perimeter.  Perhaps  as  important,  the  measured  quantities  can  refer  to  the  total  acoustic  field  rather 
than  to  its  unwanted  component  only,  and  the  methodology  can  automatically  distinguish  between  the  two. 

In  the  paper,  we  construct  a  general  solution  to  the  aforementioned  noise  control  problem.  The  appa¬ 
ratus  used  for  deriving  the  general  solution  is  closely  connected  to  the  concepts  of  generalized  potentials 
and  boundary  projections  of  Calderon’s  type.  For  a  given  total  wave  field,  the  application  of  Calderon’s 
projections  allows  us  to  definitively  split  between  its  incoming  and  outgoing  components  with  respect  to  a 
particular  domain  of  interest,  which  may  have  arbitrary  shape.  Then,  the  controls  are  designed  so  that  they 
suppress  the  incoming  component  for  the  domain  to  be  shielded  or  alternatively,  the  outgoing  component 
for  the  domain,  which  is  complementary  to  the  one  to  be  shielded.  To  demonstrate  that  the  new  noise 
control  technique  is  appropriate,  we  thoroughly  work  out  a  two-dimensional  model  example  that  allows  full 
analytical  consideration.  To  conclude,  we  very  briefly  discuss  the  numerical  (finite-difference)  framework  for 
active  noise  control  that  has,  in  fact,  already  been  worked  out,  as  well  as  some  forthcoming  extensions  of  the 
current  work:  Optimization  of  the  solution  according  to  different  criteria  that  would  fit  different  practical 
requirements,  applicability  of  the  new  technique  to  quasi- stationary  problems,  and  active  shielding  in  the 
case  of  the  broad-band  spectra  of  disturbances. 
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1.  Introduction. 


1.1.  Background.  The  area  of  active  shielding  and  noise  control  has  an  extensive  history  for  a  variety 
of  applications.  No  adequate  review  of  the  field  can  be  provided  in  the  framework  of  a  focused  research 
publication,  and  the  examples  we  mention  here  are  simply  representative.  They  do,  however,  establish  the 
relevance  of  our  study.  Elliot,  Stothers,  and  Nelson  in  [1]  develop  a  procedure  for  minimizing  the  noise  level 
at  a  number  of  pointwise  locations  (where  the  sensors  —  microphones  —  are  actually  placed);  Wright  and 
Vuksanovic  in  [2,3]  discuss  directional  noise  cancellation.  Various  issues  associated  with  the  development  and 
implementation  of  noise  cancellation  techniques  for  aircraft  industry  applications  are  analyzed  in  collection 
volumes  [4,5].  Van  der  Auweraer,  et  al.  in  [6]  studied  the  methods  of  aircraft  interior  noise  control  using 
both  acoustical  and  structural  excitation. 

The  work  by  Kincaid,  et.  al.  in  active  noise  control  [7-11]  has  been  done  under  the  assumptions  that 
the  noise  source  is  a  well  understood  monopole  with  a  few  distinct  frequencies  (in  the  200-500  Hz  range) 
and  that  the  volume  properties  of  the  supporting  medium  are  relatively  well  understood.  The  goal  was  to 
cancel  out  the  static  noise  at  a  given  collection  of  sensors  distributed  throughout  the  aircraft  cabin  by  means 
of  designing  an  appropriate  system  of  actuators.  In  references  [12]  and  [13],  Kincaid  and  Berger  solved  the 
actuator  placement  problem  for  the  mathematically  analogous  case  of  a  large  flexible  space  structure.  The 
goal  was  to  place  actuators  so  as  to  be  able  to  effectively  damp  structural  vibrations  for  large  number  of 
modes  (cf.  [14]). 

An  overview  of  the  recent  practice  in  active  control  of  noise  is  presented  in  work  by  Fuller  and  von  Flotow 
[15].  A  typical  noise  reduction  system  would  minimize  the  noise  measured  at  multiple  sensing  locations 
by  adaptively  tuning  the  control  filter.  The  filtered-a;  Least  Mean  Squares  (LMS)  adaptation  mechanism 
originally  introduced  by  Burgess  [16]  and  Widrow,  et.  al.  [17]  is  most  often  employed.  Another  adaptation 
mechanism  called  principal  component  LMS  algorithm  (PC-LMS)  presented  by  Cabell  and  Fuller  in  [18] 
offers  faster  convergence.  Good  results  can  be  achieved  at  the  sensor  locations,  but  other  locations  are  not 
considered  in  the  problem  formulation.  By  contrast,  the  approach  proposed  here  offers  the  exact  volumetric 
cancellation  of  noise. 

In  contrast  to  many  existing  noise  control  methodologies,  the  active  shielding  technique  that  we  are 
proposing  here  suggests  the  possibility  of  exact  uniform  volumetric  noise  cancellation;  moreover,  this  can  be 
done,  if  necessary,  through  the  use  of  only  surface  sensors  and  actuators.  This  paper  contains  a  systematic 
description  of  the  mathematical  foundations  of  the  new  methodology  in  the  continuous  formulatioUj  including 
a  model  example  for  the  two-dimensional  Helmholtz  equation.  The  corresponding  finite- difference  framework 
has  been  previously  introduced  and  studied  in  the  series  of  publications  by  Ryaben’kii  and  co-authors  [19-22]. 

Our  technique  actually  allows  one  to  split  the  total  acoustic  field  into  two  components,  the  first  deemed 
“friendly”  and  called  sound,  the  second  deemed  “adverse”  and  called  noise  (see  Sections  1.2,  3.2),  and 
subsequently  provides  for  the  exact  volumetric  cancellation  of  the  adverse  component  while  leaving  the  fidendly 
one  unaQected  The  input  data  for  the  noise  control  system  are  some  quantities  measured  on  the  perimeter 
of  the  region  to  be  shielded  (see  Section  2).  A  very  important  feature  of  our  approach  is  that  the  measured 
quantities  can  refer  to  the  total  acoustic  ffeld  rather  than  its  adverse  component  only.  Then,  the  foregoing 
split  into  sound  and  noise  and  subsequent  noise  cancellation  is  performed  by  the  control  system  automatically 
with  no  special  attention  required  to  it. 

Let  us  mention  that  the  analysis  tools  that  we  employ  for  building  our  methodology  involve  boundary 
integrals  of  the  kind  routinely  used  in  the  classical  potential  theory  (see  Section  2)  and  then  extend  to  the 
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operator  representation  and  generalized  boundary  projection  operators  of  Calderon’s  type  (see  Section  3). 
The  apparatus  of  the  classical  potential  theory  as  applied  to  the  Helmholtz  equation  has  been  used  in  the 
past  for  the  analysis  of  active  noise  control  problems,  see,  e.g.,  the  monograph  by  Nelson  and  Elliot  [23].  In 
this  book,  the  authors  construct  special  control  sources  on  the  surface  of  the  region  to  be  shielded  for  what 
they  call  active  absorption  and  reflection  of  the  undesirable  noise.  These  sources  are  obtained  as  dipole  and 
monopole  layers  (for  absorption)  and  only  monopole  layers  (for  reflection),  with  the  responses  in  the  form  of 
double  and  single  layer  potentials  and  only  single  layer  potentials,  respectively.  The  fundamental  diQerence 
between  our  study  and  that  of  [23]  is  that  the  authors  of  [23]  need  to  know  the  actual  boundary  trace  of  the 
noise  ffeld  to  be  cancelled  and  construct  their  control  sources  with  the  explicit  use  of  these  data.  In  other 
words,  their  analysis  in  this  part  reduces  to  the  standard  solution  of  the  boundary-value  problems  for  the 
Helmholtz  equation  using  boundary  integral  equations  of  the  classical  potential  theory.  As  opposed  to  [23], 
we  do  not  need  to  know  the  actual  “adverse”  component  of  the  acoustic  field  that  we  want  to  cancel;  this 
is  most  advantageous  because  this  component  is  obviously  impossible  to  measure  directly  in  the  presence 
of  other  acoustic  signals.  What  is  clearly  possible  to  measure  directly  is  the  overall  acoustic  field  that  may 
as  well  include  a  “friendly”  sound  component  (generated  by  the  interior  sources)  that  is  supposed  to  be 
unaffected  by  the  controls.  This  overall  ffeld  is  all  we  need  to  know  for  building  the  controls.  As  has  been 
mentioned,  our  control  methodology  is  capable  of  automatically  distinguishing  between  the  two  components 
(friendly  and  adverse)  and  accordingly  respond  only  to  the  adverse  part.  Besides,  we  provide  for  a  closed 
form  general  solution  for  the  control  sources  and  can  analyze  more  complex  settings  compared  to  the  original 
Helmholtz  equation,  in  particular,  anisotropies,  material  discontinuities,  and  certain  types  of  nonlinearities. 

Another  issue  to  be  emphasized  regarding  out  work  is  somewhat  counterintuitive  and  deviates  from  the 
“conventional  wisdom.”  It  turns  out,  in  fact,  that  one  needs  to  know  relatively  little  to  be  able  to  cancel 
the  undesirable  noise;  in  particular,  no  detailed  information  about  either  structure  or  location  or  strength 
of  the  noise  sources,  as  well  as  volumetric  properties  of  the  supporting  medium,  is  required.  We  also  stress 
that  our  methodology  yields  a  closed  form  of  the  general  solution  to  the  noise  control  problem;  availability  of 
this  general  solution  provides,  in  particular,  powerful  means  for  optimization.  Moreover,  the  same  technique 
will  apply  to  media  with  different  acoustic  conductivities  and  inhomogeneities  (within  the  linear  regime),  as 
both  the  adverse  noise  that  we  need  to  cancel  and  the  output  of  the  control  system  propagate  across  the 
same  medium.  Finally,  the  methodology  can  be  applied  to  quasi- stationary  regimes  and  in  the  future,  to 
broad-band  spectra  of  frequencies. 

1.2.  Formulation  of  the  problem.  Let  H  be  a  given  domain,  H  C  n  =  2  or  n  =  3  (higher 
dimensions  n  can  be  analyzed  similarly,  we,  however,  restrict  ourselves  by  the  cases  that  originate  from 
physical  applications).  The  domain  0  can  be  either  bounded  or  unbounded;  in  the  beginning  we  will  assume 
for  simplicity  that  12  is  bounded.  Let  T  be  the  boundary  of  Q:  F  =  d^},.  Both  on  H  and  on  its  complement 
we  consider  the  time-harmonic  acoustic  field  u  governed  by  the  non- homogeneous  Helmholtz  equation: 

Lu  =  Au  -f  k'^u  —  /.  (1-1) 

The  sources  /  in  equation  (1.1)  can  be  located  on  both  Q.  and  its  complement  to  emphasize  the 

distinction,  we  denote 


(1.2) 


3 


where  the  sources  /+  are  interior,  supp  /+  C  and  the  sources  /“  are  exterior,  supp  /  C  Accord- 

ingly,  the  overall  acoustic  field  u  can  be  represented  as  a  sum  of  two  components: 

u  —  u'^  u~ ,  (1-3) 

where 

Lu+  =  /+,  (l-4a) 

Lu-  =  /-.  (l-4b) 

Note,  both  u+  and  are  defined  on  the  entire  E”,  the  superscripts  “+”  and  refer  to  the  sources  that 
drive  each  of  the  field  components  rather  than  to  the  domains  of  these  components.  The  setup  described 
above  is  schematically  shown  in  Figure  1.1. 


Hereafter,  we  will  call  the  component  of  (1.3),  see  (1.4a),  sound  or  “friendly”  part  of  the  total  acoustic 
field;  the  component  u~  of  (1.3),  see  (1.4b),  will  be  called  noise,  or  “adverse”  part  of  the  total  acoustic  field. 
In  the  formulation  that  we  are  presenting,  will  be  a  (predetermined)  region  of  space  to  be  shielded.  This 
means  that  we  would  like  to  eliminate  the  noise  inside  H  while  leaving  the  sound  component  there  unaltered. 
In  the  mathematical  framework  that  we  have  adopted,  the  component  u~  of  the  total  acoustic  field,  i.e., 
the  response  to  the  adverse  sources  f~  (see  (1-2),  (1.3),  (1.4)),  will  have  to  be  cancelled  on  fl,  whereas  the 
component  i.e.,  the  response  to  the  friendly  sources  will  have  to  be  left  intact  on  H.  A  physically 
more  complex  but  conceptually  easy  to  understand  example  that  can  be  given  is  that  inside  the  passenger 
compartment  of  an  aircraft  we  would  like  to  eliminate  the  noise  coming  from  the  propulsion  system  located 
outside  the  aircraft  fuselage  while  not  interfering  with  the  ability  of  the  passengers  to  listen  to  the  inflight 
entertainment  programs  or  simply  converse. 

The  concept  of  active  shielding  (noise  cancellation)  that  we  will  be  discussing  implies  that  the  aforemen¬ 
tioned  goal  is  to  be  achieved  by  introducing  additional  sources  g,  suppg  C  so  that  the  total  acoustic 

field  ft  be  now  governed  by  the  equation 

Lu=f++r+  g,  (1.5) 
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and  coincide  with  only  sound  component  on  the  domain  ft: 


U  =  li"*"  , 

a;€n  xGJI 


(1.6) 


X  =  {xi,X2)  for  n  =  2  and  x  =  {xiyX2yXs)  for  n  =  3, 

The  new  sources  g  of  (1.5)  will  hereafter  be  referred  to  as  controls.  Let  us  note  that  in  practical  settings, 
when  one  typically  has  to  deal  with  acoustic  fields  composed  of  the  signals  with  multiple  frequencies  rather 
than  time-harmonic  single-frequency  fields,  active  strategies  of  noise  reduction  are  often  combined  with  the 
passive  ones.  Passive  strategies  that  employ  different  types  of  sound  insulation  have  proven  efficient  for 
damping  the  high-frequency  signals,  while  low-frequency  noise  components  more  easily  lend  themselves  to 
reduction  by  active  techniques.  We  also  note  that  the  words  “noise  cancellation”  should  not  be  interpreted 
incorrectly.  As  will  be  seen  from  the  forthcoming  analysis,  elimination  of  the  component  u~  inside  fl  also 
implies  changing  the  total  acoustic  field  outside  ft.  Therefore,  in  terms  of  the  acoustic  energy,  it  means 
redistribution  rather  than  cancellation. 

An  obvious  solution  to  the  foregoing  noise  control  problem  is  g  ^  -f~.  As,  however,  will  become 
clear,  this  solution  is  excessively  expensive.  One  one  hand,  this  expensiveness  relates  to  the  informational 
considerations  as  the  solution  g  =  -/”  requires  an  explicit  and  detailed  knowledge  of  the  structure  and 
location  of  the  sources  /“,  which,  in  fact,  will  prove  superfluous.  On  the  other  hand,  the  implementation 
of  this  solution  may  encounter  most  serious  difficulties.  In  the  previous  example  associated  with  an  aircraft, 
it  is  obviously  not  feasible  to  directly  counter  the  noise  sources  which  are  aircraft  propellers  or  turbofan  jet 
engines  located  on  or  underneath  the  wings.  Therefore,  other  solutions  of  the  control  problem,  besides  the 
most  obvious  one,  may  be  preferable  from  both  theoretical  and  practical  standpoints.  The  general  solution 
of  the  control  problem  formulated  in  this  section  is  constructed  in  the  forthcoming  Section  2. 

2.  General  solution  of  the  control  problem.  Let  us  first  introduce  fundamental  solutions  to  the 
Helmholtz  operator.  By  definition,  these  are  solutions  to  the  nonhomogeneous  equation  (1.1)  driven  by  the 
<5-source  located  at  the  origin.  For  n  =  2  we  have 

G{x)  = (2.1a) 


where  is  the  Hankel  function  of  the  second  kind  defined  by  means  of  the  Bessel  functions  Joiz)  and 

Yo{z)  as  follows  H^\z)  =  Jo{z)  —  iYo{z),  and  for  n  =  3 


^(a.)  =  - 


A'ir\x\ 


(2.1b) 


Let  us  note  that  both  solutions  (2.1)  satisfy  the  so-called  Sommerfeld  radiation  condition  at  infinity,  see, 
e.g.,  [24,25].  This  condition,  that  for  a  given  function  u{x)  is  formulated  as 

«(x)  =  O  -f  ifcu(x)  =  o  as  [x] — >  oo  (2.2a) 

for  two  space  dimensions  and 

u(x)  =  O  (|a:|“^)  ,  +  iku{x)  =  o  (|x|“^)  as  [x]  — >  oo  (2.2b) 

o\x\ 

for  three  space  dimensions,  specifies  the  direction  of  wave  propagation  and  essentially  distinguishes  between 
the  incoming  and  outgoing  waves  at  infinity.  The  Sommerfeld  condition  (2.2a)  or  (2.2b)  is  required  to 
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guarantee  uniqueness  of  the  solution  to  the  Helmholtz  equation;  unlike  for  the  Laplace  equation,  the  condition 
of  boundedness  (n  =  2)  or  vanishing  (n  =  3)  at  infinity  is  not  sufficient  for  the  the  solution  of  the  Helmholtz 
equation  to  be  unique.  Provided  that  the  following  convolution  exists  in  some  sense,  the  solution  u{x)  to 
the  nonhomogeneous  equation  (1.1)  defined  on  M"  and  satisfying  the  Sommerfeld  condition,  is  given  by 

=  f  -  y)^y- 

Next,  according  to  the  standard  procedure  (see,  e.g.,  [24]),  we  construct  the  surface  and  volume  potentials 
for  the  Helmholtz  operator  and  write  the  Green’s  formula 


-V - "  ' - - - ^ " 

volume  potential  two  surface  potentials 

which  holds  for  any  sufficiently  smooth  function  u{x).  All  integrals  on  the  right-hand  side  of  (2.4)  are, 
of  course,  convolutions,  and  n  is  the  outward  normal  to  the  boundary  P.  The  second  component  on  the 
right-hand  side  of  (2.4)  is  the  sum  of  two  surface  potentials  and  thus  it  satisfies  the  homogeneous  Helmholtz 
equation  on  H: 

The  sound  component  w+  of  the  total  acoustic  field  satisfies  equation  (1.4a)  on  M”,  and  as  the  Sommerfeld 
condition  is  built  into  the  structure  of  G{x)  (see  (2.1)),  the  function  u^{x)  can  be  obtained  as  follows 

u+{x)  =  /  Gf^dy  =  [  GLu^dy,  a:  €  E".  (2.6) 

Jn  Jn 

Then,  the  application  of  the  Green’s  formula  (2.4)  to  u+  of  (2.6)  immediately  yields 

Let  us  now  recall  that  according  to  (1.2)  and  (1.4b)  for  x  ^  there  will  be 


Lu 


=  0. 


(2.8) 


Thus,  applying  the  Green’s  formula  (2.4)  to  the  total  acoustic  field  u  =  w+  -t-  m  and  using  relations  (2.7) 
and  (2.8)  we  obtain 


r  ^  .  rf,^  ^9G  d{u++u-)^\, 

u{x)  =  GL{u+  +  u-)dy  +  (^(u+  +  w")^  -  = 


j0L.Uy  +  J(u-^ 

V - - /  V - 

U+ 


~  dn 

U~ 


X  G  Q. 


(2.9) 


Therefore,  if  we  define  the  annihilating  function  u(a:),  a;  e  fl,  as  follows 


^£o)<is. 


-u  (x). 


(2.10) 
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then  inside  Q  we  will  have 


Therefore,  if,  by  some  means,  we  obtain  the  two  quantities:  u  and  boundary  F,  then  we  can 

produce  the  annihilating  function  v  by  formula  (2.10)  and  cancel  the  adverse  component  u~  of  the  acoustic 

diL 

field  inside  fl.  In  practical  settings,  the  quantities  u~  and  -- —  or  similar  ones  should  be  actually  measured. 
A  remarkable  property  of  the  function  v{x)  defined  by  formula  (2.10)  is  that  by  virtue  of  (2.7) 


dG  du  \ 

u- - —G  dSy,  rr  G  12. 

on  on  J 


Thus,  the  quantities  to  be  measured  at  the  boundary  F,  u  and  — ,  may  refer  to  the  total  ffeld  u  =  u‘^  -\-u~ 
rather  than  its  adverse  component  u~  only.  The  annihilating  function  v  defined  as 


dG  du  \ 


automatically  filters  out  the  contribution  from  the  friendly  sound  and  responds  only  to  the  adverse  noise 
component  u~ . 

Formula  (2.13)  that  defines  the  annihilating  function  v{x)  contains  surface  integrals  of  the  quantities  u 

du/  ^  ^  ^ 

and  — .  We  will  now  discuss  the  alternative,  more  apparent,  means  of  building  v{x).  Consider  a  sufficiently 
Qn 

smooth  function  w{x)  defined  on  that  satisfies  the  Sommerfeld  condition  (2.2a)  or  (2.2b)  and  additional 
boundary  conditions  on  F: 


dw  du 

dn  dn  -  ’ 


here  u  is  the  actual  total  acoustic  field,  as  before.  In  particular,  w{x)  may  be  compactly  supported  near  the 
boundary  F.  Obviously  (because  of  (2.14)), 


v{x)  —  —  /  1  tt 


dG  du  ^ 
^  dn  dn 


G  ds,j  =  ~  {  w 


dG  dw 
^  dn  dn 


G  I  dsy,  X  G 


On  the  other  hand,  if  we  apply  the  Green’s  formula  (2.4)  to  the  function  w,  we  obtain 


.(x)  -  ^  GLwdy  =  dsy, 


An  important  consideration  now  is  to  remember  that  the  Sommerfeld  boundary  condition  (2.2a)  or  (2.2b) 
guarantees  the  uniqueness  and  therefore 

u;(a;)=:  [  GLwdy,  x  €  E".  (2.17) 

If  w{x)  is  compactly  supported,  then  integration  in  (2.17)  can  actually  be  performed  over  supp  w  rather  than 
entire  E”.  Relations  (2.16)  and  (2.17)  together  imply  that 

/■  GX«,Js  -  L  ^ 
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Therefore, 


v(x)  —  — 


1. 


GLwdy 


(2.19) 


where  the  control  g  in  (2.19)  is  defined  as 

g(x)  =  —Lw\  (2.20) 

provided  that  w  satisfies  (2.2a)  or  (2.2b)  and  (2.14).  A  variety  of  choices  for  w  (under  the  conditions  (2.2a) 
or  (2.2b)  and  (2.14))  implies  that  there  is  also  a  variety  of  controls  g  (see  (2.20))  that  solve  the  problem; 
this,  in  particular,  gives  room  for  optimization. 

Proposition  2.1.  Formula  (2.20)  describes  the  entire  variety  of  the  appropriate  controls  g. 

Proof  We  have  already  seen  that  formula  (2.20)  does  give  a  solution  to  the  active  shielding  problem.  It 
remains  to  demonstrate  that  any  solution  to  this  problem  can  be  represented  in  the  form  (2.20).  Consider 
a  control  function  g,  supp^  €  such  that  the  solution  u  of  equation  (1.5)  satisfies  the  Sommerfeld 

condition  (2.2a)  or  (2.2b)  and  also  equality  (1.6).  We  need  to  show  that  3w  that  meets  conditions  (2.2a)  or 
(2.2b)  and  (2.14)  and  such  that  g  can  be  calculated  by  formula  (2.20).  We  introduce  w{x),  x  e  as  the 
solution  to  the  nonhomogeneous  equation 


-Lw  =  g  - 


(2.21) 


subject  to  the  corresponding  Sommerfeld  condition,  (2.2a)  or  (2.2b).  Equality  (1.6)  for  the  solution  u  of 
equation  (1.5)  implies  that  the  function  w  that  solves  equation  Lw  =  g  subject  to  the  same  Sommerfeld 
condition  coincides  with  —u~  on  Consequently,  the  function  w  introduced  above  satisfies 

w\  =  u|  ,  (2.22) 

and  therefore,  it  satisfies  conditions  (2.14)  as  well.  Since  supp/"*"  C  H,  the  control  g  can  be  obtained  by 
formula  (2.20).  □ 

3.  Connection  to  the  theory  of  Calderon’s  potentials.  The  foregoing  split  between  and 
(see  formulae  (2.12),  (2.13))  can  be  conveniently  described  in  terms  of  the  generalized  potentials  and  boundary 
projection  operators  of  Calderon's  type.  We  refer  here  to  the  original  work  by  Calderon  [26]  followed  by  the 
paper  by  Seeley  [27],  and  then  work  by  Ryaben’kii  [28-30],  in  which  the  actual  form  of  the  operators  used, 
in  particular,  in  this  paper  was  introduced;  a  brief  account  of  Ryaben’kii’s  work  can  also  be  found  in  book 
by  Mikhlin,  et.  al.  [31]. 

3.1.  Potential  and  projection  for  the  domain  fl.  Consider  some  function  u{x)  such  that  Lu  —  0 
for  a:  €  fl.  Then,  the  Green’s  formula  (2.4)  yields: 


u{x) 


r  f  dG  J 

)  ’■  *  ■ 


(3.1) 


We  emphasize  that  the  representation  (3.1)  of  a  function  as  a  sum  of  a  double-layer  potential  with  the 

du 

density  u\^  and  a  single-layer  potential  with  the  density  is  valid  only  for  the  functions  that  solve 

the  homogeneous  equation  Lu  =  0  on  the  domain  H.  If,  however,  we  specify  two  arbitrary  functions  on  P 
and  substitute  them  into  (3.1)  as  densities  of  the  potentials,  then  the  resulting  function  will  obviously  be  a 
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solution  of  Xu  =  0  on  fl,  but  its  boundary  values,  as  well  as  boundary  values  of  its  normal  derivative,  will 
not,  generally  speaking,  coincide  with  the  original  densities  of  the  double-layer  and  single-layer  potentials, 
respectively.  A  generalized  potential  of  Calderon’s  type  with  the  vector  density  fr  =  specified  on  T 

is  defined  by  the  following  formula 

=  J  X  €  n,  (3.2) 

which  is  similar  to  (3.1)  except  that  we  do  not  require  ahead  of  time  that  and  in  (3.2)  be  the  boundary 

values  of  some  function  that  solves  Lu  =  0  on  H:  If,  on  the  other  hand,  and  are  the  boundary  values 

of  u  and  — —  for  some  solution  to  Lu  =  0  on  Q,  then  the  Green’s  formula  (3.1)  for  this  u  can  be  written  in 
dn 

a  shortened  form: 

For  any  (sufficiently  smooth)  function  v  specified  on  Q,  we  also  define  its  vector  trace  on  F: 

and  then  introduce  the  boundary  operator  Pr  as  a  combination  of  the  potential  Pq  oi  (3.2)  and  trace  TV* 
of  (3.4): 

Pr^r=  TVPafr.  (3-5) 

Clearly,  the  operator  Pr  of  (3.5)  is  a  projection,  Pp  =  Pp.  Indeed,  V^r  •  LPq^t  =  0,  x  €  O.  Therefore,  the 
Green’s  formula  (3.3)  yields  Pq^t  =  PnTV’Pfi^F}  which  immediately  implies  that  Pp  is  a  projection. 

The  key  property  of  the  operator  Pp  of  (3.5)  is  the  following: 

Proposition  3.1.  Those  and  only  those  vector-functions  ^p  that  satisfy  the  boundary  equation  with 
projection  (BEP) 

Pr^v  ~  Cr  (^•^) 

can  he  complemented  on  ft  to  a  function  u  such  that  Lu  =  0  on  H  and  TV  u  =  ^p . 

In  other  words,  the  range  ImPp  of  the  projection  operator  Pp  given  by  (3.5)  exhaustively  characterizes  all 
those  and  only  those  boundary  functions  that  admit  a  complement  to  the  domain  ft  in  the  foregoing  sense. 

Proof  Let  u  be  a  sufficiently  smooth  function  defined  on  11,  Xu  =  0.  Then  applying  the  operator  Tr  of 
(3.4)  to  the  Green’s  representation  (3.3)  for  u,  we  obtain  the  BEP  (3.6).  Conversely,  let  equality  (3.6)  hold 
for  some  ^p.  Denote  u  =  Pf^^r,  obviously  Xu  =  0  on  H.  In  addition,  equality  (3.6)  implies  that  TVu  =  ^p. 
Thus,  we  have  obtained  the  required  complement.  □ 

Note,  although  we  have  described  the  Calderon  potentials  and  projections  using  the  language  of  surface 
integrals,  the  previously  introduced  constructions  that  use  the  auxiliary  function  w,  in  particular,  formula 
(2.16),  will  apply  here  with  no  change.  We  emphasize  that  formula  (2.16)  essentially  shows  how  to  replace 
surface  integrals  by  volume  integrals  in  the  entire  derivation.  In  particular,  given  a  vector  density  ^p  = 
($o,fi),  we  can  take  a  sufficiently  smooth  function  w{x)  compactly  supported  near  P  such  that  (cf.  (2.14)) 

Trw  =  ^ri 
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and  redefine  the  potential  (3.2)  as  follows 


=  w{x)  -  /  GLwdy 

JQ 

=  I  GLwdy^ 


X  efi. 


Obviously,  definitions  (3.2)  and  (3.8)  are  equivalent  provided  that  relations  (3.7)  hold.  The  second  equality  in 
(3.8)  holds  because  of  the  uniqueness  in  the  form  (2.17).  Let  us  emphasize  that  the  potential  obviously 
does  not  depend  on  the  speciffc  choice  of  the  auxiliary  function  w{x)  as  long  as  condition  (3,7)  is  met 

3.2.  Split  between  u'^  and  u~.  Incoming  and  outgoing  waves.  We  now  return  to  formulae  (2.10) 

—  (2.13).  For  the  annihilating  function  v(x),  x  €  11,  we  have  v  =  ~u~,  x  €  fl.  Equation  (1.4b)  along  with  the 

(  _  du-\\ 

consideration  that  supp/"  C  implies  Lu~  ==  0,  x  G  Therefore,  if  we  denote  =  I  ,  — —  I  , 


then,  obviously,  ^p  satisfies  BEP  (3.6),  Pp^p  =  ^p  »  i-e-j  fp  €  ImPp,  and  relation  (2.10)  becomes 

v  =  -Pai-,  (3.9) 

Let  us  also  denote  f p  =  ( u,  -^1  ,  here  again  u  is  the  overall  acoustic  field,  as  in  Section  2.  Then,  we  can 

\  p 

rewrite  formula  (2.12)  as  follows 

z;  =  -Pn^p  =  -Pn^p,  X  G  (3.10) 

Thus,  the  annihilating  function  v{x),  x  G  can  be  obtained  as  a  generalized  Calderon^s  potential  with 
the  density  — ^p  =  —  (u,  I  ,  and  the  second  equality  of  (3.8)  again  indicates  that  v{x)  can  be  actually 


the  density  — fp  =  —  ti,  ,  and  the  second  equality  of  (3.8)  again  indicates  that  v{x)  can  be  actually 

\  ^^/lp 

calculated  through  the  use  of  the  auxiliary  function  w{x)  that  satisffes  boundary  condition  (3.7).  Moreover, 
as  the  potential  Pq^p  depends  only  on  ^p  and  not  on  the  particular  choice  of  w{x),  we,  in  fact,  obtain  the 
entire  variety  of  controls  (2.20)  that  cancel  out  the  same  acoustic  disturbance  u~{x)  on  the  domain 

Next,  applying  the  operator  Tr  of  (3.4)  to  both  sides  of  the  second  equality  in  (3.10)  and  taking  into 


account  that  Ppff 


we  obtain 


Pp^P  =  ir 


Equality  (3.11)  means  that  the  boundary  function  ^p,  which  is  the  vector  trace  (in  the  sense  of  (3.4))  of 
the  total  acoustic  field  u{x),  does  not,  generally  speaking,  satisfy  the  BEP  (3.6)  (unless  ^p  =  ^p 
/  du'^  \ 

Er  =  I  ^'*’5  —  =  0)  l^hus  cannot,  generally  speaking,  be  complemented  on  H  to  a  solution  of  the 

V  /  p 

homogeneous  equation  Lu  =  0.  The  portion  of  ^p  that  does  admit  the  interior  complement  is  ^p  .  What,  in 
fact,  happens  is  that  the  projection  Pp  selects  only  that  portion  of  ^p  that  can  be  complemented  on  and 
the  corresponding  complement  taken  with  the  minus  sign  provides  the  annihilating  function  v  according  to 
(3.9).  Equality  (3.10)  implies  that  the  generalized  potential  Pa  is  insensitive  to  any  contribution  to  Er  ihat 
does  not  belong  to  the  range  ImPp  of  the  boundary  projection  (3.5). 

The  foregoing  discussion  implies  that  the  split  of  the  total  acoustic  field  u{x)  into  u'^{x)  and  u~{x) 
rendered  by  formulae  (2.12),  (2.13)  (see  also  (3.10))  is,  in  fact,  realized  through  the  split  of  the  boundary 


trace  fp 


and  ^p  =  (  ^  , 
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The  component  in  (3.12)  is  obtained  by  applying  the  projection  Pp  to  according  to  formula  (3.11) 
and  of  (3.12)  complements  ^p  to  the  entire 

Assume  now  that  all  we  know  is  the  trace  ^p  (in  the  sense  of  (3.4))  on  the  boundary  T  =  of  some 
solution  u  (defined  on  M"")  to  the  equation  Lu  =  f  subject  to  the  corresponding  Sommerfeld  condition  (2.2a) 
or  (2.2b).  We  do  not  know  anything  about  the  sources  /  except  that  they  guarantee  existence  of  the  solution 
in  the  class  of  functions  satisfying  the  Sommerfeld  condition.  Then,  we  perform  the  split  (3.12)  and  see  that 
can  be  complemented  on  ft  to  the  solution  of  the  homogeneous  equation  Lu'~  =  0  (because  Pp^p  =  Cp  )♦ 
Thus,  this  component  should  be  interpreted  as  the  one  driven  only  by  the  exterior  sources  (with  respect 
to  ft).  The  component  in  (3.12)  cannot  be  complemented  on  H  to  the  solution  of  the  homogeneous 
equation  according  to  Proposition  3.1  as  formula  (3.11)  implies  that  Pp^p  =  0.  Therefore,  it  satisfies  a 
non-homogeneous  rather  than  homogeneous  differential  equation  on  ft  and  thus  should  be  interpreted  as 
the  component  driven  by  the  interior  sources.  It  is  easy  to  see  that  ^p  of  (3.12)  is  driven  only  by  the 
interior  sources  because  the  contribution  of  all  exterior  sources  to  ^p  is  already  taken  into  account  by  ^p  . 
Consequently,  fjt  of  (3.12)  can  be  complemented  to  the  solution  of  the  homogeneous  equation  outside  H,  i.e., 
on  the  domain  W^\Q>.  This  exterior  complement  will,  in  addition,  satisfy  the  Sommerfeld  condition  (2.2a) 
or  (2.2b). 

Let  us  note  that  there  is  also  a  straightforward  way  of  showing  that  of  (3.12),  for  which  Pp^p  =  0, 
can  be  complemented  on  W^\ft  in  the  sense  mentioned  above.  This  is  done  through  directly  constructing  the 
Calderon’s  potential  and  projection  for  the  complementary  domain  and  we  postpone  the  corresponding 

derivation  till  Section  3.3. 

Thus,  we  have  seen  that  the  function  of  (3.12)  that  belongs  to  the  range  of  the  projection  operator 
Pp,  G  ImPp  =  PrCf ,  represents  the  part  of  the  total  trace  ^p  that  is  accounted  for  by  the 

sources  outside  ft,  and  the  part  of  (3.12)  that  belongs  to  the  kernel  of  the  projection  operator  Pp, 
^P  e  Ker  Pp  accounted  for  by  the  sources  inside  ft.  In  other  words,  with  respect  to  the 

particular  domain  n,  the  components  and  ,  see  (3.12),  represent  boundary  traces  of  the  incoming  and 
outgoing  portions,  respectively,  of  the  total  acoustic  ffeld. 

In  fact,  the  entire  space  Hp  of  vector- functions  ^p  =  (^Oj^i)  defined  on  F  =  dft  can  be  split  into  the 
direct  sum  of  the  range  and  kernel  of  the  projection  operator  Pp: 

Hp=ImPp0KerPp.  (3.13) 


Equality  (3.13)  means  that  for  any  given  fp  6  Hp  there  is  always  a  unique  representation  in  the  form  (3.12). 
Again,  the  physical  interpretation  of  equality  (3.13)  is  that  any  given  acoustic  ffeld  can  be  split  into  the 
incoming  and  outgoing  components  with  respect  to  a  given  domain  H.  This  split  can  be  performed  on  the 
boundary  only  by  applying  the  Calderon^s  projection  Pp  to  the  trace  of  the  total  acoustic  ffeld. 

To  conclude  this  section,  let  us  only  emphasize  a  useful  and  non-trivial  consequence  of  formulae  (2.14) 
and  (2.20).  Formula  (2.20)  describes  the  variety  of  controls  g{x)  through  the  flexibility  in  choosing  the 
auxiliary  function  w{x),  which  should  still,  however,  satisfy  boundary  conditions  (2.14).  Equalities  (2.14) 
can  be  rewritten  as  (cf.  (3.7)) 


du 

dn 


(3.14) 


Trw  = 

Clearly,  if  we  add  some  fp  =  ^  KerPp,  to  ^ 

we  will  recover  a  diQerent  control  function g{x) ,  but  still  obtain  the  same  annihilating  function  v{x),  because 


on  the  right-hand  side  of  equality  (3.14),  then 
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v(x)  =  —Pn 


\  f  du\ 

r 

,  and  the  potential  Pn  is  insensitive  to  any  €  KerPr-  According  to 


aforementioned  alteration,  i.e.,  addition  of  ^pj  cannot j  of  course,  add  any  new  elements  to  the  overall  set  of 
control  functions  g.  This  alteration  only  provides  an  alternative  way  of  describing  (i.e.,  parameterizing)  some 
of  the  controls,  which  may,  however,  appear  helpful  for  optimization,  especially  in  the  discrete  framework. 


3.3.  Potential  and  projection  for  the  complementary  domain  E^\n.  The  Calderon’s  potential 
Pq,  see  (3.8),  (3.7),  and  projection  Pp,  see  (3.5),  are  introduced  for  the  interior  domain  O.  In  a  very  similar 
way,  one  can  introduce  the  Calderon’s  potential  and  projection  for  the  exterior  domain  E^\n.  Given  a 
vector  density  fp  =  (^o,6)  we  again  take  a  sufficiently  smooth  auxiliary  function  w{x)  compactly  supported 
near  F  and  such  that  it  satisfies  boundary  conditions  (3.7),  and  define  the  exterior  potential  QiRn\n^p(x), 
X  G  E^\n,  as  follows  (cf.  (3.8)): 


QR"\a^r(3^)  “ '^(^)  —  /  GLwdy 

=  [  GLwdy,  xeR^\Ct. 

Jq 


(3.15) 


Again,  similarly  to  (3.8)  the  second  equality  in  (3.15)  holds  due  to  the  uniqueness  in  the  sense  of  (2.17),  and 
the  potential  does  not  depend  on  the  particular  choice  of  w{x)  as  long  as  condition  (3.7)  holds. 

The  corresponding  underlying  construction  that  uses  the  Green’s  formula  for  obtaining  the  potential 
of  (3.15)  will  be  almost  identical  to  the  one  that  we  have  used  for  obtaining  the  interior  potential  Pn,  see 
(3.2),  except  that  the  Green’s  formula  similar  to  (3.1)  will  have  to  be  written  for  the  exterior  domain  E^\n 
as  well.  The  exterior  projection  operator  Qp,  =  Qp,  is  obtained  by  applying  the  trace  Tr  given  in  (3.4) 
to  the  potential  QiRn\n  given  in  (3.15): 


Qr^r  =  Tr 


(3.16) 


Similarly  to  Proposition  3.1,  the  range  ImQp  of  the  projection  operator  Qp,  see  (3.16),  contains  those 
and  only  those  functions  $p  that  can  be  complemented  on  the  exterior  domain  E”\n  to  a  solution  of  the 
homogeneous  equation  Xu  =  0  that  would  also  satisfy  the  corresponding  Sommerfeld  boundary  condition 
(2.2a)  or  (2.2b).  Comparing  the  definitions  (3.8)  and  (3.5)  of  the  potential  Pn  and  projection  Pp  for 
the  domain  Q  with  the  definitions  (3.15)  and  (3.16)  of  the  potential  and  projection  Qr  for  the 

complementary  domain  E^\n,  we  immediately  see  that  Qp  =  I  -  Pp,  and  consequently  ImPp  =  KerQp 
and  ImQp  —  KerPp.  Therefore,  any  boundary  function  from  KerPp  can  indeed  be  complemented  of  the 
exterior  domain  E^\n  to  the  solution  of  the  homogeneous  differential  equation  (as  has  been  mentioned  in 
Section  3.2). 

Similarly  to  the  split  (3.13)  we  can  also  obtain 

Hr  =  ImQp©KerQp.  (3.17) 


The  functions  from  ImQr  andKerQp  are  boundary  traces  (in  the  sense  of  (3.4))  of  the  incoming  and  outgoing 
wave  ffelds,  respectively,  where  the  terms  “incoming”  and  “outgoing”  refer  now  to  the  exterior  domainMT-\Q. 
rather  that  interior  domain  H  as  before. 

It  is  important  to  note  that  when  constructing  both  Pq  and  Pp,  see  Section  3.2,  and  and  Qp, 

see  above,  we  use  the  same  operator  Tr  of  (3.4),  i.e.,  the  same  direction  n  of  differentiation  at  the  boundary 
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r.  This  basically  creates  an  asymmetry  because  for  one  domain  the  normal  is  directed  inwards  and  for 
another  one  —  outwards.  However,  this  asymmetry  allows  us  to  always  consider  the  same  space  of  boundary 
functions  Ep  and  thus  arrive  at  convenient  relations:  Pp  =  /  —  Qp,  and  decompositions  (3.13)  and  (3.17). 
Hereafter,  we  will  always  assume  that  the  direction  of  normal  differentiation  at  the  boundary  F  is  the  same 
for  both  the  interior  and  exterior  domains. 

Representation  (3.17)  suggests  that  one  can  also  describe  the  controls  g{x)  for  the  interior  domain  Q 
in  terms  of  the  potential  and  projection  Qp  constructed  for  the  exterior  domain  Indeed, 

KerQp  contains  those  and  only  those  boundary  functions  that  can  be  interpreted  as  traces  of  the  outgoing 
waves  with  respect  to  the  domain  M”\n.  The  sole  task  of  the  controls  g{x)  is  to  eliminate  the  influence 

that  domain  M^\n  exerts  on  the  domain  H,  or  in  other  words,  eliminate  the  component  of  the  acoustic  ffeld 

(  Q'lx  \ 

that  is  outgoing  with  respect  to  Thus,  given  the  total  field  ^p  =  fw,  j  boundary,  the 

control  has  to  respond  to  its  outgoing  component  (/  —  Qr)  ^p  ^  KerQp  (with  respect  to  the  domain  E”\n) 
and  actually  eliminate  it.  In  other  words,  it  has  to  reconstruct  —  ( J  —  Qr)ir  from  the  given  data  ^p;  the 
application  of  the  operator  -(/  -  Qp)  via  auxiliary  function  w{x)  immediately  yields  the  set  of  controls  in 
the  form  (2.20). 

We  finally  note  that  the  construction  of  Qur.\a  and  Qp,  which  is  fully  parallel  to  Pa  and  Pp,  essentially 
lifts  the  assumption  made  in  the  beginning  of  Section  1.2  that  the  domain  to  be  shielded  should  be  bounded. 

4.  A  more  general  formulation  of  the  active  shielding  problem.  Here  we  discuss  several  impor¬ 
tant  generalizations  to  the  formulation  of  the  noise  control  problem  analyzed  previously. 

4.1.  Changing  the  far- field  boundary  conditions.  Instead  of  the  entire  space  E”,  consider  a  suf¬ 
ficiently  large  domain  Qo,  H  C  Ho  £  E””.  We  introduce  the  setup  similar  to  that  of  Section  1.2.  Namely,  the 
total  acoustic  field  u{x)  is  governed  by  the  Helmholtz  equation  (1.1),  the  sources  f{x)  that  drive  the  solution 
are  again  split  into  the  interior  and  exterior  components,  see  (1.2);  supp/"*”  C  Q  as  before,  but  as  opposed 
to  Section  1.2,  supp/“  C  (rather  than  E^\r^).  Also  different  is  the  far-held  boundary  condition. 

Previously  we  required  that  the  solution  u{x)  satisfy  the  Sommerfeld  condition  (2.2a)  or  (2.2b);  now  we 
specify  a  general  linear  homogeneous  boundary  condition  at  which  we  formulate  as  an  inclusion: 

u  6  C/o,  (4.1) 

and  require  that  equation  (1.1)  subject  to  boundary  condition  (4.1)  be  uniquely  solvable  on  the  domain 
Ho  for  any  (sufficiently  smooth)  right-hand  side  f{x).  Similarly  to  Section  1.2,  we  decompose  the  overall 
solution  into  the  friendly  component  u'^{x),  x  e  fio,  which  is  called  sound,  and  adverse  component  u“(x), 
X  G  flo?  which  is  called  noise,  see  formulae  (1.3),  (1.4).  Again  similarly  to  Section  1.2,  we  would  like  to 
construct  a  control  function  g{x)y  supp^  C  so  that  the  total  acoustic  held  u{x),  which  then  includes 

the  output  of  controls  as  well  and  is  governed  by  equation  (1.5),  coincide  on  the  domain  Q  with  only  sound 
component  u'^{x),  i.e.,  satisfy  equality  (1.6).  We  denote  by  G  the  inverse  operator  to  L  so  that  u  =  G  f, 
u  gUq]  the  aforementioned  unique  solvability  means,  in  particular,  that  for  any  w  ^Uq  we  have 

w=GLw.  (4.2) 

In  the  previously  analyzed  case,  the  space  Uoy  see  (4.1),  contained  all  those  and  only  those  functions  that 
were  dehned  on  E’^  and  satished  the  corresponding  Sommerfeld  condition  (2.2a)  or  (2.2b),  and  the  inverse 
operator  G  was  constructed  as  a  convolution  with  the  fundamental  solution  (2.1a)  or  (2.1b).  As  for  the 
uniqueness  in  the  form  (2.17),  it  has  now  transformed  into  (4.2). 
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Using  the  inverse  operator  G,  we  will  now  formally  construct  the  potentials  and  projections  for  both 
domains,  and  similarly  to  how  it  has  been  done  in  Section  3.  Given  a  boundary  density  =  (?o,  6), 

we  take  some  sufficiently  smooth  auxiliary  function  w  ^  Uq  that  satisfies  condition  (3.7)  and  define 


Pa^r{x)  =  w-G{Lw\J 


(4.3a) 


Pr^r  =  TrPa(r> 


(4.3b) 


Definitions  (4.3a)  and  (4.3b)  are  obviously  analogous  to  (3.8)  and  (3.5),  respectively.  Expressions  in  curly 
brackets  in  (4.3a)  mean  that  after  applying  the  differential  operator  L  to  the  function  w{x)  the  result  is 
truncated  to  the  corresponding  domain,  either  Q,  or  Q>q\Q..  The  second  equality  in  (4.3a)  is  guaranteed  by 
the  uniqueness  (4.2). 

Similarly,  for  the  complementary  domain  Qo\^  we  have: 


Qno\(iir{x)  =  w-  G 

=  G{Lw\^},  xeno\Cl, 


(4.4a) 


QrCr  =  T'l'  QfioVn^r-  (4.4b) 

Again,  the  definitions  (4.4a)  and  (4.4b)  are  analogous  to  (3.15)  and  (3.16),  respectively,  and  the  second 
equality  in  (4.4a)  follows  from  (4.2). 

The  definitions  (4.3)  and  (4.4)  are  so  far  only  formal.  To  make  them  meaningful,  we  need  to  ensure 
that  the  projections  (4.3b)  and  (4.4b)  indeed  split  the  overall  wave  field  into  the  incoming  and  outgoing 
components,  as  before.  In  other  words,  we  need  to  guarantee  the  key  result  similar  to  that  of  Proposition  3.1. 
This,  in  fact,  requires  that  the  potentials  (4.3a)  and  (4.3b)  depend  only  on  the  density  fr  and  be  independent 
of  the  particular  choice  of  iu(a:)  as  long  as  it  satisfies  (3.7).  In  other  words,  we  need  to  require  that  if  =  0, 
then  for  any  function  w(x),  w  €  Uq,  with  the  zero  trace,  Trw  =  Cr  =  0  (see  (3.4)),  we  would  obtain: 

Pn^r  =  Pr20r  =  rc-G{lu;|j^}  =  0,  x€n,  (4.5a) 

and 

Qno\n^r  =  Qno\n0r  =  ■n' -  <3  =  0,  x  €  (4.5b) 

One  can  show  that  for  the  Helmholtz  operator  L  of  (1.1)  and  boundary  condition  (4.1)  that  guarantees 
the  unique  solvability,  equalities  (4.5)  will  always  hold.  Alternatively,  we  can  say  that  a  particular  choice 
of  the  trace  according  to  (3.4)  is  proper  for  the  given  differential  operator  L  in  the  sense  that  zero  trace 
implies  zero  potential.  The  same  is,  in  fact,  true  for  many  other  linear  elliptic  differential  operators  of 
the  second  order.  For  detail,  we  refer  the  reader  to  [29],  where  requirements  of  the  type  (4.5)  are  used 
for  characterizing  a  special  class  of  boundary  traces,  the  so-called  clear  traces,  which  are  then  employed 
for  building  the  generalized  potentials  and  projections  of  Calderon’s  type.  In  the  current  paper  we  restrict 
ourselves  to  considering  only  the  traces  of  type  (3.4)  (Cauchy  data),  which  are  applicable,  in  particular,  to 
the  second-order  time-harmonic  wave  problems  that  we  are  solving. 
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Having  constructed  the  potentials  and  projections  (4.3)  and  (4,4)  based  on  the  operator  G  that  corre¬ 
sponds  to  the  general  boundary  condition  (4.1),  we  can  now  develop  the  full  set  of  arguments  similar  to  that 
of  Section  3  provided  that  equalities  (4.5)  hold.  In  particular,  the  annihilating  function  t;(x)  for  the  interior 
domain  Q  is  still  given  by  formula  (3.10)  with  the  only  difference  that  the  potential  Pq  is  now  defined  by 
(4.3a).  There  is  still  a  split  (3.13)  of  the  total  acoustic  field  into  the  incoming  and  outgoing  components 
with  respect  to  the  domain  however,  the  split  is  now  rendered  by  the  projection  Pr  defined  by  (4.3b)  via 
the  operator  G  and  thus  differs  from  that  of  Section  3.  Similarly,  there  is  a  split  (3.17)  of  the  total  acoustic 
field  into  the  incoming  and  outgoing  components  with  respect  to  the  domain  although  the  projection 

Qr  defined  by  (4.4b)  via  the  operator  G  is  again  different  from  that  of  Section  3. 

Since  the  projection  operators  Pr  and  Qr  =  I~Pt  now  not  the  same  as  those  of  Section  3,  so  are  the 
results  Prfr  and  Qr^r  of  applying  these  operators  to  a  given  ^r-  In  other  words,  the  incoming  and  outgoing 
parts  of  a  given  $r  obtained  by  applying  the  operator  Pr  of  (4.3b)  (or  Qr  of  (4.4b))  will  notj  generally 
speaking,  be  the  same  as  the  incoming  and  outgoing  parts,  respectively,  obtained  in  Section  3.  This  means 
that  the  decomposition  of  the  total  acoustic  field  into  the  incoming  and  outgoing  components  depends  on 
the  external  boundary  condition  (4-^)  that  these  components  satisfy,  which  is,  of  course,  reasonable  from  the 
standpoint  of  physics.  However,  one  can  show  that  although  the  change  of  the  domain  and  external  far-field 
boundary  condition  introduced  here  compared  to  Section  3  does  amount  to  the  change  of  the  projection 
operators  Pr  and  Qr,  it  does  not  induce  any  changes  of  the  subspace  ImPr  =  Ker  Qr-  As  concerns  the 
other  subspace  of  the  split  (3.13)  or  (3.17),  KerPr  =  Im  Qr,  it  does,  generally  speaking,  change.  This  can 
be  interpreted  as  follows.  Im  Pr  contains  all  those  and  only  those  boundary  functions  from  ^r  that  can  be 
complemented  on  the  interior  domain  O  to  the  solution  of  the  homogeneous  equation  Xu  =  0.  This  set  of 
functions  as  a  whole  will  not,  of  course,  depend  on  any  external  boundary  condition.  However,  we  actually 
construct  this  set  by  applying  the  projection  Pr  to  the  entire  space  5r.  The  direction  of  projecting  depends 
on  the  external  boundary  conditions.  This  dependency  manifests  itself  through  the  change  in  the  operators, 
and  for  a  given  ^r  €  Sr  we  will  generally  speaking  obtain  a  different  PrCr,  which  will  still  belong,  however, 
to  the  same  subspace  Im  Pr  that  does  not  change.  The  change  in  the  direction  of  projecting  onto  this 
subspace,  ImPr,  obviously  implies  the  change  in  the  second  component  of  the  direct  sum  (3.13)  (or  (3.17)). 
This  second  component  Ker  Pr  =  Im  Qr  contains  all  those  and  only  those  functions  from  Hr  that  can  be 
complemented  on  r2o\f^  to  a  solution  of  Lu  =  0  that  satisfies  (4.1)5  these  functions  will,  of  course,  depend 
on  the  actual  exterior  domain  and  far-filed  boundary  condition. 

Let  us  now  emphasize  the  most  important  thing  regarding  the  relation  between  the  potential  and  pro¬ 
jection  operators  that  we  build  according  to  formulae  (4.3)  and  (4.4)  (cf.  (3.2),  (3.5)  and  (3.15),  (3.16))  and 
the  noise  control  problem  that  we  analyze  using  these  operators.  It  turns  out  that  although  the  potentials 
and  projections  change  with  the  change  of  the  extetnal  boundai^  condition  and  so  does  the  decomposition 
of  the  acoustic  ffeld  into  the  incoming  and  outgoing  components  (with  respect  to  or  Qo\Q>),  the  set  of 
control  functions  {g{x)}  that  we  obtain  is  essentially  insensitive  to  these  changes.  Indeed,  as  the  annihilat¬ 
ing  function  u(x)  is  given  by  (3.10)  with  the  new  potential  defined  by  (4.3a),  then  similarly  to  (2.20)  and 
Proposition  2.1  the  entire  set  of  appropriate  controls  for  Q.  is  given  by 


g{x)  -Lw 


(4.6) 


where  w{x)  satisfies  boundary  condition  (4.1)  and  also  (2.14)  (the  latter  can  be  rewritten  as  (3.14)).  Thus, 
the  only  difference  between  (2.20)  and  (4.6)  is  that  in  (2.20)  w{x)  is  supposed  to  satisfy  one  of  the  Sommerfeld 
conditions  (2.2a)  or  (2.2b)  rather  than  condition  (4.1).  On  the  other  hand,  for  applications  we  will  primarily 
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be  interested  in  controls  concentrated  near  the  perimeter  of  the  region  to  be  shielded,  i.e.,  near  the  boundary 
r  =  rather  than  far  away  from  it.  Therefore,  we  can  always  consider  only  those  w{x)  that  are  compactly 
supported  near  F,  Such  functions  will  obviously  satisfy  any  homogeneous  far-field  boundary  condition  and 
thus  the  set  of  controls  built  on  the  basis  of  these  compactly  supported  tt;(a:;)’s  will  be  exactly  the  same  if  we 
use  either  of  the  formulae  (2.20)  or  (4.6). 

4.2.  Changing  the  differential  operator.  Instead  of  the  original  Helmholtz  operator,  see  (1.1),  we 
will  now  consider  some  other  linear  differential  operator  that  operates  on  the  functions  defined  on  the  domain 
Ho.  We  will  keep  the  same  notation  L  for  simplicity  and  require  that  the  acoustic  field  u{x)  governed  by  the 
new  equation  Lu  =  /  satisfy  the  same  far- field  boundary  condition  (4.1)  as  before.  Similarly  to  Section  4.1, 
we  will  require  that  the  equation  Lu  —  f  subject  to  boundary  condition  (4.1)  be  uniquely  solvable  for  any 
sufficiently  smooth  right-hand  side  /,  and  also  for  simplicity,  we  will  keep  the  same  notation  G  for  the 
corresponding  inverse  operator.  From  the  standpoint  of  the  physical  model,  the  change  of  the  operator  may 
be  accounted  for,  e.g.,  by  the  change  in  the  properties  of  the  medium,  across  which  the  acoustic  signals 
propagate.  For  example,  as  opposed  to  the  previously  analyzed  case  of  the  classical  constant-coefficient 
Helmholtz  equation  that  describes  stationary  waves  in  an  isotropic  medium,  we  may  consider  anisotropies  of 
the  medium,  e.g.,  obstacles  like  passenger  seats  and  overhead  bins  that  can  be  introduced  in  the  passenger 
compartment  of  an  aircraft.  Basically,  obstacles,  i.e.,  anisotropies,  of  various  nature  can  be  introduced  in 
either  of  the  domains  ft  or  or  in  both  of  them.  At  any  rate,  provided  that  the  regime  of  the  wave 

propagation  is  still  linear  (some  questions  related  to  nonlinearities  are  touched  upon  in  Section  4.4)  we  arrive 
at  a  new  equation  that  may  have  variable  coefficients  in  the  areas  where  the  obstacles  are  introduced.  In 
this  section,  we  limit  our  analysis  by  the  case  when  these  coefficients  are  still  smooth  across  the  boundary 
F;  the  case  when  discontinuities  in  the  material  properties  are  allowed  along  F  is  considered  in  Section  4.3. 

Using  the  new  operators  L  and  G  we  can  repeat  the  entire  analysis  of  Section  4.1.  In  so  doing,  we 
arrive  at  the  operators  Pn,  Pr  and  Qr  defined  by  the  same  formulae  (4.3)  and  (4.4),  respectively, 

but  based  on  the  new  L  and  G.  Provided  that  relations  (4.5)  hold,  we  can  guarantee  the  results  similar 
to  those  of  Proposition  3.1  and  thus  decompose  the  overall  acoustic  field  into  the  incoming  and  outgoing 
components.  (The  terms  ‘incoming”  and  ‘^outgoing”  in  this  section  still  relate  to  a  particular  domain,  H 
or  but  the  filed  itself  and  thus,  each  of  its  components,  is  governed  by  a  different  equation.)  As  has 

been  mentioned  (see  Section  4.1  and  also  [29]  for  details),  relations  (4.5)  essentially  imply  that  the  trace 
operator  Tr  is  correlated  in  some  sense  with  the  differential  operator  L.  In  this  paper  we  are  using  only 
the  traces  of  type  (3.4),  which  are  consistent  (see  [29]),  in  particular,  with  the  second-order  linear  elliptic 
differential  operators  T,  for  which  the  Neumann  boundary  data  reduce  to  a  standard  normal  derivative,  i.e., 
Lu  =  V  (pVit)  +  {lower  order  terms},  p  =  p{x)  is  a  scalar  function.  Thus,  the  changes  of  the  operator  that 
we  are  considering  shall  only  be  within  this  class.  This  is  obviously  going  to  be  sufficient  for  all  practical 
purposes  related  to  time- harmonic  acoustics. 

Of  course,  the  potentials  and  projections  built  on  the  basis  of  the  new  operators  L  and  G  will  not  be 
the  same  as  those  of  either  Section  4.1  or  Section  3.  Similarly  to  the  considerations  of  Section  4.1,  one  can 
see  that  changing  L  only  in  the  exterior  domain  fto\ft  will  change  both  Pr  and  Qr  ImPr  =  KerQr  will 
remain  the  same,  whereas  KerPr  -  ImQr  will,  generally  speaking,  change.  (The  change  of  KerPr  =  ImQr 
accounts,  as  before,  for  a  different  direction  of  projecting  onto  ImPr.)  Changing  L  only  on  the  interior 
domain  will  again  cause  changes  to  both  Pr  and  Qr  but  ImQr  =  KerPr  will  not  change  (whereas 
KerQr  =  ImPr  will).  However,  neither  of  these  changes  will  essentially  aQect  the  control  functions  g(x). 
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Indeed,  to  cancel  out  the  unwanted  noise  on  the  interior  domain  we  use  the  annihilating  function 
v(x)  =  -Pn^r  with  (310))  and  accordingly,  arrive  at  the  control  functions  g(x) 

given  by  formula  (4.6),  where  L  is  now  the  new  operator  and  is  a  sufficiently  smooth  function  that 
satisfies  (4.1)  and  (2.14).  Formula  (4.6)  still  gives  the  general  solution  for  controls  and  formally  looks  the 
same  as  before  with  the  “hidden”  difference  in  plugging  in  the  new  L.  The  reason  that  the  controls  are 
described  by  similar  formulae  even  when  the  actual  wave  propagation  processes  are  governed  by  different 
equations  is  that  in  every  particular  instance  both  the  unwanted  disturbances  to  be  cancelled  and  the  output 
of  controls  propagate  across  the  same  medium  and  satisfy  the  same  far-ffeld  boundary  conditions.  It  turns 
out  that  having  realized  it  we  can  design  the  controls  with  no  explicit  knowledge  of  the  actual  properties  of 
the  supporting  medium,  Le,,  the  coeflcients  of  L  except  in  the  region  where  w{x)  ^  0,  i.e.,  where  we  apply 
the  operator  X,  see  (4.6). 

As  has  been  mentioned,  the  region  where  w{x)  ^  0  may  always  be  chosen  as  a  narrow  strip  straddling 
the  boundary  F.  (This  is  beneficial  from  the  standpoint  of  practical  design.)  Therefore,  in  the  particular 
case  when  the  aforementioned  obstacles  (either  in  Q.  or  in  flQ\Q.  or  in  both  domains)  are  introduced  not 
right  next  to  the  boundary  but  rather  at  a  distance,  the  operator  L  near  the  boundary  remains  the  same 
Helmholtz  operator  of  (1.1),  and  thus  the  exact  same  control  functions  as  those  obtained  in  Sections  2  and 
3  on  the  basis  of  compactly  supported  if;(a;)’s  can  be  used  for  noise  cancellation  in  the  current  more  general 
setting. 

4.3.  Introducing  discontinuities  in  the  material  properties.  The  next  generalization  that  comes 
naturally  after  considering  the  changes  of  the  far-field  boundary  conditions  (Section  4.1)  and  the  changes  of 
the  coefficients  of  L  that  may  be  caused  by  anisotropies  in  Q.  and/or  no\n  (Section  4.2),  is  to  consider  two 
different  media  separated  by  the  boundary  F.  In  some  sense,  this  is  the  ultimate  change  of  the  operator  that 
we  may  consider;  it  further  extends  the  analysis  of  Section  4.2,  in  which  we  allowed  for  the  variation  of  the 
coefficients  of  L  but  still  assumed  that  they  were  smooth  across  the  boundary  F.  In  this  section,  we  allow 
the  coefficients  of  the  operator  L  to  have  jumps  on  the  interface  F  (e.g.,  the  media  on  different  sides  of  F 
may  have  different  refraction  indices  k,  see  (1.1)).  Thus,  some  additional  interface  conditions  may,  generally 
speaking,  be  required  for  solvability.  In  many  cases,  and  in  particular  those  that  we  are  studying  in  this 
paper,  these  additional  conditions  are  the  continuity  of  the  solution  itself  and  the  continuity  of  its  normal 
flux  across  the  interface.  For  the  second-order  operators,  we  can  write  it  in  the  general  form  as  follows: 

^  (4.7) 

where  and  are  traces  of  the  solution  u{x)  in  the  sense  of  (3.4)  taken 

from  the  interior  and  exterior  sides  of  F,  respectively,  and  is  a  linear  operator  that  provides  the  required 
relation  between  these  traces  (e.g.,  continuity  of  the  solution  and  its  normal  flux  on  F). 

We  require  that  the  differential  equation  Lu  =  f  be  uniquely  solvable  on  the  domain  Q>q  for  any  suffi¬ 
ciently  smooth  right-hand  side  /  in  the  class  of  functions  that  satisfy  the  far-field  boundary  condition  (4.1) 
and  the  interface  compatibility  condition  (4.7).  We  denote  by  G  the  corresponding  inverse  operator  to  L 
so  that  u  =  Gf,  ueUo  and  u  satisfies  (4.7).  The  aforementioned  unique  solvability  implies,  in  particular, 
that  relation  (4.2)  will  hold  for  any  function  w{x)  that  satisfies  (4.1)  and  (4.7). 

To  build  the  Calderon’s  operators,  we  now  need  to  distinguish  between  the  boundary  traces  on  different 
sides  of  F.  Namely,  given  the  boundary  density  we  take  an  arbitrary  (sufficiently  smooth)  auxiliary 
function  w{x)  €  Uq,  so  that  it  also  meets  the  interface  condition  (4.7)  (in  the  sense  that  the 
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two  functions  =  Tri^o\n)^  connected  via  (4.7))  and  construct  the  operators  Pq  and 

Pr  according  to  formulae  (4.3);  the  boundary  trace  in  formula  {4.3b)  is  taken  from  the  side  of  the  domain 
Q,.  Again,  the  condition  of  (4.5a)  type  is  required  to  guarantee  the  result  of  Proposition  3.1.  Similarly, 
we  obtain  Qqq\q  and  Qr  according  to  formulae  (4.4);  these  operators  operate  on  and  the  auxiliary 

function  w{x)  €  Uq,  needed  for  the  calculation  is  again  taken  so  that  it  meets  the 

interface  condition  (4.7).  Proposition  3.1  will  hold  provided  equality  (4.5b)  is  satisfied. 

One  important  difference  compared  to  the  previously  analyzed  cases  is  that  here,  generally  speaking, 
I  —  Pp  Qp  and  therefore,  we  can  no  longer  claim  that  ImPp  =  KerQp  and  KerPp  =  ImQp.  In  fact, 
similarly  to  what  we  had  before  all  those  and  only  those  boundary  functions  fp^^  that  satisfy  the  BEP 
i.e.,  G  ImPp,  can  be  complemented  on  H  to  the  solution  of  the  homogeneous  equation 
Lu  =  0  (Proposition  3.1);  these  functions  G  ImPp  should  still  be  interpreted  as  traces  of  the  incoming 
waves  with  respect  to  the  domain  Q.  Moreover,  the  projection  operator  Pp  still  renders  the  split  of  the 
overall  acoustic  field  into  the  direct  sum  of  its  incoming  and  outgoing  components  with  respect  to  the 
domain  fl.  However,  contrary  to  what  we  had  before  we  now  need  to  treat  this  overall  acoustic  field  ^p^^ 
as  defined,  i.e.,  “measured,”  on  the  interface  P  from  its  interior  side  and  thus,  rewrite  formula  (3.13)  as 
=  ImPp  0  KerPp.  Moreover,  the  functions  G  KerPp  that  are  interpreted  as  traces  of  the  outgoing 
waves  with  respect  to  H,  cannot,  generally  speaking,  be  complemented  to  the  solution  of  the  homogeneous 
equation  Xu  =  0  on  HoXH,  i.e.,  do  not  belong  to  ImQr-  This  “mismatch”  between  KerPp  and  ImQp  may 
give  rise  to  reflections  from  the  interface  when  analyzing  a  particular  problem. 

Of  course,  the  considerations  that  we  have  just  brought  forward  are  fully  reciprocal  with  respect  to 
the  exterior  domain  In  other  words,  we  rewrite  the  decomposition  formula  (3.17)  as  = 

ImQp  0  KerQp  emphasizing  that  we  now  measure  the  overall  solution  from  the  exterior  side  of  the  interface 
r,  and  note  that  the  outgoing  waves  ^  KerQp  may  be  (partially)  reflected  from  the  interface  because 

KerQp  ImPp. 

As  concerns  the  set  of  controls  {g{x)}  however,  it  still  remains  essentially  the  same  as  before  and  neither 
of  the  aforementioned  changes  in  the  formulation  of  the  problem  and  solution  structure,  including  possible 
reflections  from  the  interface,  actually  affect  it.  Indeed,  the  annihilating  function  for  the  domain  0.  is  given 

du^ 


by  v{x)  =  -Pn$r^\  where  = 


u, 


dn 


(H) 


,  see  (3.10),  i.e.,  we  use  the  generalized  Calderon’s  potential 


with  the  density  obtained  as  a  trace  of  the  overall  acoustic  field  measured  on  the  interior  side  of  F.  Using 
equalities  (4.3a),  we  rewrite  it  as 


(4.8) 


where  the  auxiliary  function  w{x)  belongs  to  Uq  see  (4.1),  =  fp^\  and  also  satisfies  the  interface 

condition  (4.7).  Accordingly,  the  control  functions  g{x)  is  again  given  by  formula  (4.6),  where  we  can  always 
take  w{x)  compactly  supported  near  F.  Moreover,  from  representation  (4.8)  for  the  annihilating  function 
we  conclude  that  we  do  not  even  need  to  know  the  entire  w{x)y  it  is  sufficient  to  construct  this  auxiliary 
function  only  on  the  exterior  domain  making  sure  that  i.e.,  that 

it  satisfy  the  correct  interface  compatibility  condition  (4.7).  In  practice,  the  trace  of  the  solution  at  the 
boundary  represents  the  values  that  are  actually  measured.  Therefore,  we  are  building  the  controls  according 
to  formula  H.6),  where  the  auxiliary  function  w{x)  is  Huned”  in  the  sense  of  to  the 

trace  of  the  actual  acoustic  ffeld  measured  on  the  exterior  side  ofV. 

Moreover,  representation  (4.8)  for  the  annihilating  function  v{x)  along  with  formulae  (4.4)  and  (3.17) 
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suggest  that  this  annihilating  function  is,  in  fact,  constructed  so  that  to  eliminate  the  outgoing  waves  with 
respect  to  the  domain  Q,q\Q..  This  interpretation  has  already  been  mentioned  in  the  end  of  Section  3.3, 
however  then  we  analyzed  the  case  when  the  outgoing  waves  with  respect  to  r2o\n  were  the  same  as  the 
incoming  waves  with  respect  to  (7.  In  the  current  framework  this  is  no  longer  so  and  the  waves  generated 
by  the  sources  /~  not  only  enter  the  domain  ft  but  may  also  be  reflected  back  from  the  interface  F.  The 
controls  that  we  build  according  to  formula  (4.6)  respond  to  the  actual  filed  measured  on  the  exterior  side  of 

cance/  out  the  entire  outgoing  jf eld  component  with  respect  to  Qo\^;  both 
the  portion  of  the  ffled  that  propagates  through  the  boundary  F  into  Q,  and  the  portion  that  may  be  reflected 
back  to  fto\Q.  Note,  the  reflections  themselves  do  not  “contaminate”  the  measurements  because  they  already 
belong  to  ImQr  and  thus  do  not  affect  the  output  of  the  controls.  We  conclude  that  in  some  sense  it  may 
be  more  natural  to  describe  the  controls  in  terms  of  the  outgoing  waves  with  respect  to  rather  than 

incoming  waves  with  respect  to  0,  especially  in  the  case  when  these  waves  are  not  the  same. 

4.4.  Introducing  nonlinearities.  A  recipe  for  including  nonlinearities  in  the  formulation  of  the  active 
shielding  problem  that  we  analyzed  in  this  paper  was  first  proposed  by  Ryaben’kii  in  [22].  In  this  section, 
we  simply  mention  with  no  detail  that  the  last  argument  of  Section  4.3  actually  suggests  that  nonlinearities 
of  a  particular  nature  can  be  handled  by  the  same  controls.  Indeed,  the  controls  are  designed  so  that  to 
cancel  out  the  entire  outgoing  component  of  the  acoustic  field  with  respect  to  a  particular  domain.  Clearly, 
as  long  as  the  nonlinearities  are  not  concentrated  on  or  immediately  near  the  perimeter  of  the  region  to  be 
shielded,  the  controls  (4.6)  will  still  be  responding  to  the  actual  acoustic  field  measured  at  the  boundary  and 
the  nature  of  the  sources  and  the  supporting  medium  away  from  the  boundary  —  whether  they  are  linear  or 
nonlinear  —  will  not  matter.  This  conclusion  holds  provided  that  the  nonlinearities  satisfy  one  important 
limitation,  namely,  that  they  do  not  alter  the  frequency  of  the  acoustic  signals.  It  is  known,  however,  that 
many  nonlinearities  do  produce  multiple  frequencies  from  a  single-frequency  input.  In  this  case,  provided 
that  the  spectrum  of  frequencies  is  available  ahead  of  time,  the  controls  can  still  be  built  as  described  above, 
but  for  each  frequency  separately  and  independently. 

5.  An  example.  The  synthesis  of  our  control  input  will  be  demonstrated  on  a  simple  two-dimensional 
example,  where  the  circular  region  Cl  defined  by  r  <  R  is  protected  from  exterior  noise  by  active  control 
acting  either  along  the  perimeter  or  within  an  annular  region  surrounding  Cl.  We  derive  the  exact  analytic 
solution  in  Fourier  representation,  and  give  explicit  expressions  for  the  optimal  distributed  control.  No 
information  about  the  nature  of  the  exterior  sound  sources  is  required,  since  the  solution  depends  only  on 
the  sound  field  and  its  normal  derivative  measured  along  the  perimeter.  In  actual  implementation,  this  exact 
control  function  would  be  approximated  by  a  finite  number  of  acoustic  inputs.  The  same  approach  applies 
to  complicated  geometries,  which  can  be  handled  numerically. 

5.1.  Helmholtz  equation  in  E^.  The  Helmholtz  equation  (1.1)  for  n  =  2  expressed  in  polar  coordi¬ 
nates  {r,0)  reads 

Urr  +  Ur/r  +  uee/r"^  +  k’^u  =  /,  (5.1) 

where  u  is  the  amplitude  of  the  time- harmonic  acoustic  field  u  e^^^,  and  k  =  c/co  is  the  wavenumber  (refraction 
index). 

5.2.  Fourier  transform.  Let  u  be  the  Fourier  transform 

27r 

ue-^”^^d9 


19 


and 


1  ^27r 

The  Helmholtz  equation  (5.1)  in  this  Fourier  representation  becomes 

Urr  -I-  Ur/v  +  (/c^  -  VT? It^)u  =  /. 


(5.2) 


The  functions  u  and  /  are  associated  with  a  particular  mode  number  m.  This  is  implicit  in  our  notation  so 
that  the  equations  are  not  cluttered  with  an  excessive  number  of  superscripts. 

5.3.  Boundary  conditions.  Functions  u  must  be  analytic  at  the  origin  r  =  0  and  decay  to  zero  at 
infinity.  Moreover,  we  shall  consider  compactly  supported  forcings  /  so  that  as  r  — ^  oo,  the  responses  u 
must  represent  outgoing  waves.  The  outgoing  waves  satisfy  the  Sohimerfeld  condition  (2.2a),  which  means 
that  asymptotically 

u(r)  ^  as  r  — ^  oo  (5.3) 

y/r 

for  some  constant  a. 

5.4.  Fundamental  solution.  Let  us  compute  the  response  to  a  unit  monopole  input  at  r  =  s,  ^  =  0, 
i.e.,  a  shifted  fundamental  solution.  (This  presents  no  loss  of  generality,  because  for  a  particular  location  of 
the  source  we  can  always  achieve  that  its  angular  coordinate  6  be  zero  by  rotating  the  coordinate  system.) 
We  have  a  ^-type  source 


f{r^0)  =  6{r  —  s)S{rO)  =  (5(r  —  s)6{0)/r  =  S{r  —  s)S{0)/s 


so  that 

=  (5^4) 

The  solution  of  equation  (5.2)  with  the  right-hand  side  (5.4)  is  explicitly  constructed  as  follows.  We  first 
notice  that  equation  (5.2)  is  homogeneous  everywhere  except  at  r  =  s.  The  homogeneous  counterpart 
to  equation  (5.2)  has  two  linearly  independent  eigensolutions,  which  we  take  as  =  Jm{kr)  and 

=  Y{kr)->riJm{kr)  =  iH^'>{kr),  where  Jm(kr)  and  Ym{kr)  are  the  Bessel  functions,  so  that  w^‘^(r) 
satisfies  the  condition  of  analyticity  at  r  =  0  mentioned  in  Section  5.3  and  vP\r)  satisfies  the  condition 
(5.3).  Then,  we  build  the  solution  of  the  nonhomogeneous  equation  (5.2)  driven  by  the  5-source  (5.4)  from 
the  two  branches;  (r)  for  r  •<  s  and  (r)  for  r  ^  s,  so  that  in  addition  it  be  continuous  at  r  =  s 
and  have  a  particular  discontinuity  in  the  first  derivative  prescribed  by  (5.4).  In  so  doing,  we  arrive  at  the 
fundamental  solution 


G(r,  s) 


1 

(^^) (-^^^ (^^)  "b  ^*^m(^^))  V  S 

];Jm{ks){Ym{kr)  +  iJm{kr))  for  r  >  s 
V  4 


(5.5) 


5.5.  Subspaces  of  incoming  and  outgoing  waves.  Let  us  now  consider  a  distribution  of  sources 
on  the  plane  E^.  Then,  in  the  Fourier  space  we  obtain  equation  (5.2)  driven  by  the  right-hand  side  f{r), 
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which  has  an  extended  support  on  r  >  0,  as  opposed  to  the  point-wise  support  like  in  (5.4).  In  this  case, 
the  corresponding  solution  is  given  by  the  integral  (27rs  is  the  Jacobian): 

u(r)=  j  G{r,s)f{s)27tsds  ,  (5.6) 

supp/ 

which  is  obviously  a  “linear  combination”  of  the  fundamental  solutions  G{r,  s)  centered  at  different  locations 
on  supp  /,  and  which  can  be  considered  as  a  generalization  of  the  standard  convolution  with  the  fundamental 
solution  (cf.  (2.3))  that  we  have  in  the  case  of  constant  coefficients. 

At  a  given  location  R,  we  can  now  easily  distinguish  between  the  two  components  of  the  overall  solution 
u{r)  given  by  (5.6):  u‘^(r)  generated  by  the  interior  sources  /  (s),  s  <  and  u~{r)  generated  by  the 
exterior  sources  /  (s)^  s  >  R.  Indeed,  similarly  to  Section  2  we  have 

ti+(r)  =  J  G{r,s)f'^{s)27rsds  =  J  G{r,s)f{s)27Tsds  ,  (5.7a) 

r+  {s<R} 

supp/ 


u  (r)  = 


G{r^s)f  {s)27rsds' 


G{r,  s)f{s)  27rs  ds 


Therefore,  we  can  say  that  those  and  only  those  solutions  u‘^(r)  that  can  be  attributed  to  the  interior 
sources  in  the  sense  of  representation  (5.7a)  are  “parallel”  to  the  right  branch  w^^^(r)  =  Y{kr) +iJm{kr)  of 
the  fundamental  solution  (5.4),  i.e.,  the  following  Wronskian  is  equal  to  zero: 


du^ 

.  dr 

dr  . 

Analogously,  those  and  only  those  solutions  u  (r)  that  can  be  attributed  to  the  exterior  sources  (with  respect 
to  the  location  r)  in  the  sense  of  representation  (5.7b)  are  “parallel”  to  the  left  branch  u^^\r)  =  Jm{kr)  of 
the  fundamental  solution  (5.4),  i.e.,  the  following  Wronskian  is  equal  to  zero: 


u~ 

du~ 

dw(i)  “  ^  • 

(5.8b) 

-  dr 

dr  J  r-R 

Clearly,  the  portion  u'^{r)  of  the  total  solution  u{r)  which  is  entirely  due  to  the  interior  sources  s  <  R 
satisfies  the  homogeneous  counterpart  of  the  differential  equation  (5.2)  on  the  exterior  domain  t  >  R  and 
thus,  using  a  standard  property  of  the  Wronskian  that  it  is  identically  zero  once  it  is  zero  at  one  point,  we 
obtain  from  (5.8a): 

<^0^rnikr) +iJmikrj)^+  _  ^  =0,  r>R.  (5.9a) 

Similarly,  the  portion  u~{r)  due  to  the  exterior  sources  s>  R  alone  satisfies 


dJ-fn(kr) 

- i - Lii 

dr 


^  .  U/U/  „  _ 

u  -Jmikr)-^  =  0,  r  <  R. 
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Note,  relations  of  type  (5.8)  and  (5.9)  have  been  used  extensively  for  constructing  the  so-called  artificial 
boundary  conditions  that  are  needed  for  the  numerical  solution  of  infinite-domain  problems,  see  the  review 
paper  by  Tsynkov  [32]. 

The  kernels  of  the  two  linear  constraints  (5.9a)  and  (5.9b)  define  a  decomposition  of  the  solution  space 
into  a  direct  sum  of  subspaces  corresponding  to  outgoing  and  incoming  waves,  respectively.  Note,  the  concept 
of  splitting  the  solution  into  incoming  and  outgoing  components  was  introduced  and  explained  in  detail  in 
Section  3.2.  Once  we  identify  the  location  r  =  R  with  the  boundary  of  the  domain  we  can  use  relations 
(5.9)  to  actually  decompose  of  the  overall  acoustic  field  u{r)  into  the  interior  and  exterior  contributions, 

u^{r)  and  u”(r),  with  respect  to  H.  This  is  done  by  solving  equations  (5.9)  along  with  u'^  +  u~  =  u  and 

“t"  d'LL 

- 1 - =  — ,  i.e.,  a  total  of  four  linear  equations,  with  respect  to  the  unknown  quantities 

dr  dr  dr 


while  treating  u  and  —  as  the  given  data.  The  solution  is 
dr 

=  T  (yrniM+iJmikr)) 


■KT  [  du  ^  ^  ^dJm{kr)  \  d{Ym{kr)  +  iJm{kr)) 

dr  J  dr 


u^\ 

\  1 

4- 

dr  1^ 

U  1 

1  Trr 

\r=R  ~  2 

du~ 

7rr 

dr 

r=R  2 

-  ^iY^{kr)+iJ,n{kr))\Jm{kr) 


-  -riYm{kr)  +  iJm{kr)) 
dr 


dJm 

dr 


(5.10a) 


(5.10b) 


(5.10c) 


(5.10d) 


5.6.  Noise  control  along  the  perimeter.  In  this  section,  we  build  the  controls  concentrated  only  on 
the  boundary  of  the  domain  H,  i.e.,  on  the  circle  r  —  R,  We  consider  the  situation  where  the  total  acoustic 
field  is  generated  by  three  categories  of  sources:  interior  sources  /  ,  monopole-type  control  inputs  g  along 
the  perimeter  r  =  R,  or  in  other  words,  as  we  are  discussing  the  Fourier  representation,  control  function 
p(r)  A5{r  —  R)  with  the  point- wise  support  at  r  =  i?,  and  exterior  sources  /  .  The  responses  to  these 
inputs  will  be  v  and  fi”,  respectively.  The  governing  equation  is  (1-5).  In  Fourier  representation,  the 
total  response  is  denoted  by  u  =  u'^  +  v u~ .  Note,  in  the  general  discussion  of  Section  2  we  have  also  first 
obtained  the  annihilating  function  v{x)  as  a  response  to  only  surface  excitation,  see  formulae  (2.10),  (2.12), 
(2.13),  and  then  recovered  the  same  function  as  an  output  of  the  volumetric  control  (2.20)  constructed  with 
the  help  of  the  auxiliary  function  w{x).  In  the  framework  of  the  current  example,  we  analyze  volumetric 
controls  in  Section  5.9. 

The  interior  region  can  be  completely  shielded  from  the  exterior  noise  by  devising  control  inputs  which 
exactly  compensate  for  the  exterior  contribution  to  it,  such  that  equation  (2.11)  holds.  This  condition  is 
equivalent  to  requiring  that 

u~{R)  +  v{R)  =  0  (5.11) 


/  du  dv\ 

hm  ^ 

rT-R  \  dr  dr  J 
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Both  V  and  u~'  are  responses  to  exterior  sources  relative  to  the  open  domain  so  both  v  and  u~  separately 
satisfy  equation  (5.9b).  Since  (5.9b)  and  (5.11)  imply  (5.12),  the  condition  (5.12)  is  redundant.  Therefore, 
(5.11)  is  a  sufficient  condition  which  we  will  use.  Integral  (5.6)  along  with  the  representation  g{r)  =  AS{r-R) 
for  the  control,  in  which  the  constant  A  is  yet  to  be  determined,  immediately  yields  that  for  r  <  R: 
v{r)  =  A  •  27rR  •  G{r,R),  and  using  the  continuity  of  v{r)  and  equality  (5.11),  we  conclude  that  A  = 
-u~  {R)/{27rR  •  G{R,  R)).  Consequently,  the  control  input  is  given  by  the  formula 


..  . _ ^u"{R) _ Hr-R) 

Jm{kR){Ym{kR)+iJmikR))  2-nR 


(5.13a) 


which  can  be  evaluated  using  equation  (5.10c)  and  simplified  to  read: 


^  d\og{Ym{kr)  +  iJmjkr)) 
\  dr  dr 


5{r  -  R) 

r=R 


(5,13b) 


In  compliance  with  the  general  theory  developed  in  Sections  2  and  3,  expression  (5.13b)  requires  the  knowl¬ 
edge  of  the  total  u  and  du/dr  along  the  exterior  of  the  perimeter  r  =  R,  but  does  not  require  any  information 
about  the  nature  of  exterior  sound  sources.  Moreover,  the  control  (5.13b)  exactly  shields  the  interior  region 
r  <  R  from  the  exterior  noise,  while  leaving  the  interior  sound  component  completely  unaffected.  Realizable 
approximations  of  this  exact  solution  may  prove  an  effective  method  of  reducing  exterior  noise. 


5.7.  Control  of  a  single  exterior  source.  As  has  been  shown  in  Sections  2  and  3,  and  also  follows 

from  the  previous  analysis  in  Section  5.5,  the  control  g  given  by  (5.13b)  is  insensitive  to  the  contribution 

du 

of  interior  sources  to  u  and  — .  Let  us  consider  a  single  exterior  source  located  aX  s  >  R  whose  forcing 

dr 

function  /  is  given  by  (5.4).  For  this  exterior  source  the  response  function  is  u~{R)  =  (5(R,  s),  see  (5.6), 
and  substituting  this  expression  into  (5.13a),  we  obtain 


__  Ymjks)  -\-iJm{ks)  S{r  -  R) 
Ym{kR)  +  iJm{kR)  2nR 


(5.14) 


Exterior  noise  field  before  control 


Knowing  the  asymptotic  behavior  of  Bessel 

functions  as  m  — >>  oo,  we  can  immediately  see 

,  .  .  fR\'"S{r-R)  , 

that  g  =  g^r^[  —  ]  _  „  and  therefore 

\s  J  ZttH 

when  s  >  R  the  Fourier  series  with  the  coeffi¬ 
cients  g^  is  convergent  to  an  infinitely  smooth 
function  of  the  argument  0, 

5.8.  A  computational  excunple.  As  a 
computational  example,  we  have  considered 
the  problem  of  shielding  a  domain  of  radius 
R  =  1  from  the  time- harmonic  acoustic  dis¬ 
turbances  with  the  wavenumber  k  —  5  (see 
(1.1)).  A  particular  exterior  noise  field  u~  was 
generated  by  the  forcing  function 


Fig.  5.1.  A  particular  sound  ff eld  u  (r,  0)  generated  by  the  exte¬ 
rior  sources  given  by  formula  (5.15).  The  vertical  scale  is  in  decibels. 


f  {x,y)  =  S{x  -  3)5{y  -  1) 

(5.15) 

-S{x  -  5)^(2/  +  1) 

which  produces  the  noise  field  u~  depicted  in 
Figure  5.1. 
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1^1  Perimeter  control  g 


Fig.  5.3.  The  sound  ffeld  v{r,9)  generated  by  the  control  input 
along  the  perimeter.  The  vertical  scale  is  in  decibels. 


Noise  field  with  control 


Fig.  5.4.  The  sound  ffeld  u~{r,6)  +  v{r,9)  resulting  after  the 
perimeter  control  has  been  applied.  The  vertical  scale  is  in  decibels. 


Figure  5.2  shows  the  required  control  in¬ 
put  along  the  perimeter,  which  can  be  com¬ 
puted  using  equation  (5.13b). 

This  perimeter  input  produces  a  control 
acoustic  field  v  depicted  in  Figure  5.3. 

With  control  applied,  the  noise  level 
within  the  protected  region  r  <  i?  =  1  is  re¬ 
duced  to  zero.  This  can  be  seen  in  Figure  5.4. 

We  note  that  in  Figures  5.1,  5.3,  and 
5.4  of  this  section  the  vertical  scale  is  given 
in  decibels;  these  units  are  defined  as  db  = 
201ogio  1^1,  where  A  is  the  amplitude  of  a  sig¬ 
nal,  or  alternatively,  db  =  lOlogio  where 
\A\‘^  is  power. 

5.9.  Optimal  control  with  annular 
support.  As  opposed  to  Section  5.6,  in  which 
we  analyzed  surface  controls,  here  we  con¬ 
sider  volumetric  control  sources  compactly 
supported  on  the  annular  region  of  thickness  a. 
In  other  words,  we  now  have  interior  sources 
in  the  region  r  <  R,  exterior  sources  in  r  >  J^, 
and  control  sources  within  R  <  r  <  R  a. 
Our  goal  will  be  to  find  a  control  which  can¬ 
cels  out  the  unwanted  exterior  noise  on  and 
is  also  optimal  in  the  sense  of  a  cost  criterion 
specified  below  (see  (5.17)). 

We  note  that  according  to  Proposition  2.1 
the  entire  variety  of  appropriate  controls  is  de¬ 
scribed  by  formula  (2.20),  where  the  auxiliary 
function  w{x)  satisfies  the  Sommerfeld  con¬ 
dition  (2.2a)  and  boundary  conditions  (2.14). 
To  obtain  compactly  supported  controls  it  is 
natural  (as  has  been  repeatedly  pointed  out 
in  Sections  2,  3,  and  4)  to  consider  compactly 
supported  functions  w{x).  However,  this  ap¬ 
proach  will  not,  generally  speaking,  yield  all 
possible  compactly  supported  controls  because 
a  function  w{x),  which  is  not  compactly  sup¬ 
ported  near  the  boundary,  may  nonetheless 
generate  a  compactly  supported  g{x)  accord¬ 
ing  to  formula  (2.20).  Still,  we  expect  that  for 
the  problems  that  originate  from  applications 
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this  approach  will  on  one  hand  provide  for  the  only  tangible  way  to  parameterize  the  compactly  supported 
controls  and  on  the  other  hand,  the  set  of  controls  obtained  this  way  will  be  sufficiently  wide  for  subsequent 
optimization. 

For  the  purpose  of  the  demonstration  in  this  section  it  will  be  convenient  to  consider  a  different  sub-class 
of  compactly  supported  controls,  namely  square  integrable  controls.  The  simple  formulation  that  we  study 
here,  i.e.,  the  formulation  that  presumes  the  isotropy  of  the  supporting  medium  everywhere  and  provides  for 
the  separation  of  variables  along  the  boundary,  will  allow  us  to  explicitly  construct  the  optimal  control  in 
this  class,  where  the  optimality  is  treated  as  a  minimum  of  the  standard  L2  norm,  see  (5.17). 

Using  the  same  argument  as  before,  we  see  that  the  interior  will  be  exactly  shielded  from  the  exterior 
sources  provided  that  the  linear  constraint  (5.11)  holds  at  the  boundary  r  —  R.  This  constraint  written 
explicitly  is 


u  (R)  -h 


G{R,  s)g{s)  27vsds  —  0 


which  is  satisfied  by  many  alternative  distributions  of  control  sources.  Subject  to  this  linear  constraint,  we 
shall  seek  to  minimize  the  cost  function  defined  by  the  L2  norm  of  the  control  effort: 


pR-[-a 

llsf  =  /  |5(s)|2  27rsds 

JR 


Let  us  consider  the  space  of  square  integrable  control  inputs  g  \  [R,  R a]  C.  At  the  optimal  solution  g, 
which  is  a  point  in  this  space,  the  corresponding  level  set  of  the  control  effort  norm  must  be  tangent  to  the 
hyperplane  of  complex  codimension  1  defined  by  the  linear  constraint  (5.16).  If  we  perturb  the  optimal  g  by 
a  perturbation  7  within  the  constraint  set,  the  cost  must  never  decrease: 


f 

JR 


(Iff (s)  p  +  l7(s) p  +  9{s)j{s)  +  ff  (s)7(s))  2wsds>\\g\{- 


This  can  only  happen  if  all  allowable  perturbations  7  satisfy  the  following  (real)  orthogonality  condition: 


-I  pR-^-Of  ____ 

(ff,  7)  2  (.9{sh(s)  +  g{s)'y{s))  2ns ds  =  0 


Since  the  perturbation  7  is  limited  only  by  the  single  complex  constraint  (obtained  as  a  variation  of  (5.16)): 

/R-\-(i 

G{R,  5)7(5)  2ns  ds  =  0  ,  (5.19) 


the  orthogonality  condition  (5.18)  imposes  a  severe  limit  on  the  form  of  g.  We  note  that  (5.19)  implies  that 
all  allowed  perturbations  7(5)  are  orthogonal  to  all  functions  of  the  type  p{s)  =  C  •  (j(it!,  s),  where  (7  is  a 
complex  constant,  in  the  sense  of  (5.18): 

1  pR-^a, 

(p,j)  =  ^  /  {CGiR,  s)7(s)  +  CG{R,  3)7(5))  2nsds  =  0,  VC  G  C  ,  (5.20) 

2  Jr 

because  the  integral  in  (5.20)  is  the  real  part  of  the  integral  in  (5.19)  premultiplied  by  C.  Since  C  is 
arbitrary,  the  converse  also  holds,  so  that  conditions  (5.19)  and  (5.20)  are  equivalent.  Therefore,  the  complex 
one-dimensional  subspace  of  functions  {p{r)}  =  {C  •  G{R,r)\C  G  C}  spanned  by  the  function  G{R,r)  is 
orthogonal  to  all  allowed  perturbations.  It  follows  from  (5.18)  that  the  optimal  control  g  must  lie  within  the 
subspace  {p].  Therefore,  within  the  annular  support  of  the  control  input  R<r<R^-a,  g{r)  =  C  •  G{R,  r) 
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for  some  complex  constant  C  (elsewhere,  g{r)  =  0).  As  the  last  step,  the  constant  of  proportionality  C  is 
obtained  from  (5.16).  The  unique  solution  for  the  control  that  we  obtain  by  this  method  is  given  by 


5(0  =  — 


u-{R)  G{R,r) 

pR-\-a, 

/  \G{R,  s)\'^  27r sds 

Jr 


4u  {R)  {Y^jkr)  -  iJmjkr)) 


(5.21) 


7TJm{^R)^s‘^[Jm{ksy  -  Jm~l{ks)Jm+l{ks)  +Ym{ks)‘^  -  Ym~l{ks)Ym+l{ks)]^ 


s=R-\-a 


for  <  r  <  +  a,  and  g{r)  =  0  otherwise.  The  term  u~{R)  in  (5.21)  must  be  evaluated  via  the  equation 

(5.10c)  which  requires  the  measurement  of  u  and  du/dr  along  the  perimeter  r  =  R  prior  to  the  application 
of  the  control.  Alternatively,  this  process  of  measurement  and  control  can  be  iterated  in  order  to  adapt  to 
slow  changes  in  the  noise  field.  The  minimum  cost  of  this  optimal  control  is  \\g\\^  =  \u~{R)\^. 

Returning  to  the  computational  example 
of  Section  5.8,  the  optimal  distributed  com 
trol  input  to  exterior  noise  sources  given  by 
formula  (5.15)  was  computed  using  equation 
(5.21).  In  this  example  we  used  R  —  1  and 
a  =  3.  The  amplitude  \g\  of  the  resulting  op¬ 
timal  control  is  depicted  in  Figure  5.5. 

6.  Discussion  and  conclusions.  We 
have  presented  an  accurate  mathematical  for¬ 
mulation  of  the  problem  of  active  shielding  of 
a  predetermined  region  of  space  from  time- 
harmonic  acoustic  disturbances.  We  have  con¬ 
structed  a  general  solution  of  this  problem 
in  the  closed  form,  and  using  the  appara¬ 
tus  of  generalized  Calderon’s  potentials  and 
boundary  projection  operators  analyzed  sev¬ 
eral  consecutively  more  complex  cases:  “Pure” 
Helmholtz’  equation  on  the  entire  space  with 
the  Sommerfeld  boundary  conditions  at  infinity,  other  types  of  the  homogeneous  far- field  boundary  conditions 
that  may  be  set  either  at  infinity  or  at  a  finite  external  boundary,  spatial  anisotropies  and  discontinuities  in 
the  material  properties,  and  certain  kinds  of  nonlinearities. 

Our  approach  to  the  problem  of  active  noise  control  possesses  several  key  advantages.  It  does  not 
require  any  knowledge  of  either  structure  or  location  or  strength  of  the  actual  sources  of  noise  that  is 
about  to  be  cancelled.  Neither  does  it  require  knowledge  of  the  properties  of  the  medium  across  which 
the  acoustic  signals  propagate,  except,  maybe  right  next  to  the  boundary  of  the  domain  to  be  shielded. 
It  guarantees  the  exact  volumetric  cancellation  of  the  unwanted  noise  throughout  the  domain  of  interest, 
while  the  input  information  sufficient  for  building  the  controls  is  given  by  a  particular  set  of  measurements 
performed  only  at  the  perimeter  of  this  domain.  Moreover,  the  control  sources  themselves  are  concentrated 
also  only  on  or  near  the  perimeter  of  the  region  to  be  shielded  thus  rendering  an  effective  surface  control 
of  the  volumetric  properties.  Finally,  the  controls  are  constructed  so  that  to  eliminate  only  the  component 


Optimal  distributed  control  input 


Fig,  5.5.  The  magnitude  of  the  optimal  distributed  control  input 
|g(r,  0)1  with  support  1  <  r  <  4  for  the  computational  example  of 
Section  5.8  with  noise  sources  given  by  formula  (5.15). 
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of  the  acoustic  field  generated  by  the  exterior  sources  (i.e.,  located  outside  the  domain  of  interest),  while 
leaving  the  component  that  comes  from  the  interior  sources  completely  unaffected.  In  so  doing,  however, 
the  measurements  performed  on  the  perimeter  of  the  domain  can  refer  to  the  total  acoustic  field  rather  than 
its  unwanted  exterior  component  only,  and  the  methodology  can  automatically  tell  between  the  interior  and 
exterior  contributions  to  the  overall  field. 

To  demonstrate  that  the  technique  is  appropriate,  we  have  thoroughly  analyzed  a  model  two-dimensional 
example,  and  using  the  separation  of  variables  in  the  Helmholtz  equation  written  in  polar  coordinates  and 
the  apparatus  of  Bessel  functions,  constructed  both  purely  surface  and  near-surface  volumetric  controls  for 
shielding  a  cylindrical  region  from  a  given  distribution  of  the  noise  sources. 

In  the  future,  the  discrete  framework  for  the  active  noise  control  is  going  to  be  used  for  analyzing 
complex  configurations  that  originate  from  practical  designs.  This  discrete  framework  has,  in  fact,  already 
been  developed  to  a  large  extent,  see  [19-22].  It  possesses  the  same  set  of  attractive  features  as  the  foregoing 
continuous  formulation.  The  discrete  formulation  is  based  on  the  difference  potentials  method  by  Ryaben’kii 
[28—30],  and  uses  finite-difference  analogues  of  Calderon’s  potentials  and  boundary  projection  for  constructing 
the  controls  and  analyzing  their  properties.  As  opposed  to  the  continuous  model  described  in  this  paper,  in 
the  finite-difference  framework  both  the  measurements  are  performed  and  the  control  sources  are  located  at 
the  grid  nodes,  i.e.,  there  are  discrete  sets  of  sensors  and  actuators  for  the  acoustic  field,  which  is  obviously 
more  feasible  from  the  standpoints  of  physics  and  engineering  applications. 

There  are  most  important  issues  yet  to  be  addressed,  in  particular  in  the  discrete  framework,  aimed  at 
eventually  creating  practical  designs.  First,  there  are  considerations  of  conditioning  —  specifically,  how  the 
measurement  errors  that  are  inevitable  in  any  practically  engineered  system  of  sensors  will  propagate  through 
the  control  system.  Second,  there  are  considerations  of  optimality.  As  we  have  seen,  there  is  an  entire  family 
of  control  functions  that  can  eliminate  the  unwanted  noise  on  a  given  domain  —  this  is,  in  fact,  true  for 
both  continuous  and  discrete  formulations.  For  a  simple  model  example  in  the  continuous  formulation,  we 
have  shown  in  Section  5.9  how  to  find  a  particular  representative  of  this  family,  which  is  optimal  in  the  sense 
of  a  certain  criterion.  Generally,  optimizing  the  general  solution  for  controls  is  a  separate  substantial  task 
composed  of  a  number  of  sub-problems. 

Namely,  the  criteria  for  optimization  (i.e.,  objective  functions)  that  would  fit  different  practical  require¬ 
ments  need  to  be  clearly  identified.  These  criteria  will  certainly  be  problem-dependent.  For  example,  the 
designer  of  a  noise  suppression  system  in  the  passenger  compartment  of  an  aircraft  should  obviously  try 
and  minimize  its  weight  and  energy  consumption,  while  for  suppressing  the  household  appliances’  noise  the 
primary  concern  may  be  the  cost  of  the  active  control  system.  After  identifying  the  objective (s)  for  opti¬ 
mization,  an  appropriate  optimization  methodology  has  to  be  chosen;  it  may  be  either  gradient-based  or 
combinatorial  (or  a  combination  of  the  two)  and  may  also  include  the  initial  effectiveness  filters.  Moreover, 
there  may  be  different  levels  for  optimization.  In  the  beginning,  we  can  search  through  the  family  of  exact 
solutions  to  the  noise  control  problem  that  is  at  our  disposal  (see  Sections  2,  3,  and  4)  and  thus  obtain  the 
one,  which  will  be  optimal  in  some  sense.  In  so  doing  we  are  still  guaranteed  that  whatever  optimum  we  find 
it  will  still  do  the  job,  or  in  other  words,  exactly  cancel  out  the  unwanted  noise  on  the  domain  of  interest. 
If,  however,  the  optimum  we  find  this  way  is  still  unsatisfactory  from  the  standpoint  of  the  chosen  criteria, 
we  can  go  beyond  the  optimization  of  the  exact  solution  only.  In  fact,  we  can  introduce  the  tolerance,  i.e., 
acceptable  level  of  noise  reduction  as  opposed  to  the  exact  cancellation,  and  using  this  greater  flexibility  try 
and  better  satisfy  the  optimization  criteria.  This  is  already  an  approximate  optimization  because  it  deviates 
from  the  exact  solution  for  controls;  as  such,  it  will  certainly  require  solving  the  finite-difference  governing 
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equation(s)  repeatedly  inside  the  optimization  loop  for  the  purpose  of  assessing  the  quality  of  the  noise 
reduction.  (Again,  in  the  previous  case  the  noise  cancellation  is  exact  and  this  is  known  ahead  of  time.)  A 
prerequisite  of  this  optimization,  of  course,  is  an  efficient  analysis  code;  most  likely,  fast  parallel  codes  will 
be  needed  for  solving  real-life  problems. 

In  the  paper,  we  have  addressed  only  the  case  of  time-harmonic  disturbances.  In  fact,  the  methodology 
will  work  for  quasi-stationary  acoustic  fields  as  well,  i.e,,  those  that  change  slowly  in  time  (slowly  on  the  scale 
of  the  inverse  frequency  of  time-harmonic  oscillations) .  In  the  latter  case,  practical  implementation  of  the 
noise  suppression  technique  will  require  real-time  measurements.  If,  however,  we  are  dealing  with  stationary 
interior  and  exterior  fields,  we  can  use  the  simplified  version  of  the  technique  that  permits  precalculation  of 
the  dependence  of  the  response  on  domain  geometry  and  system  properties.  In  other  words,  a  microphone 
is  needed  in  the  final  system  only  to  sense  phase  and  amplitude.  Together  with  this,  preprogrammed  noise¬ 
cancelling  actuators  are  sufficient. 

Finally,  we  should  mention  that  time-harmonic  (i.e.,  Helmholtz)  problems  of  active  shielding  (as  well  as 
the  aforementioned  quasi-stationary  problems)  represent  only  a  portion  of  the  overall  variety  of  formulations 
that  the  potential  users  of  such  methodologies  would  like  to  have  explored.  Even  though  these  Helmholtz 
problems  are  of  a  substantial  significance  themselves  from  the  standpoints  of  both  mathematics  and  ap¬ 
plications,  it  is  clear  that  true  time- dependent  problems,  i.e.,  those  that  accommodate  broadband  acoustic 
fields,  will  ultimately  be  demanded  by  practitioners  and  therefore  deserve  theoretical  and  numerical  study. 
The  problem  of  active  noise  control  in  the  formulation  that  involves  broadband  spectra  of  frequencies  was 
studied  by  Malyuzhinets  [33]  and  Fedoryuk  [34].  Some  initial  results  on  the  time-dependent  active  shielding 
problem  obtained  in  the  framework  of  generalized  Calderon’s  potentials  can  be  found  in  work  by  Zinoviev 
and  Ryaben’kii  [35].  Generally,  the  extension  to  broadband  spectra  of  disturbances  is  nontrivial  in  every 
respect  and  provides  a  novel  challenge  from  the  standpoints  of  both  theoretical  analysis  and  subsequent 
practical  implementation. 
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